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INTRODUCTION 


Background 

Current  airfield  pavement  design  procedures  do  not 
include  a  methodology  to  determine  the  thickness  of  a 
stabilized  base  course  material.  In  October  of  198 ,  a 
Federal  Aviation  Administration  (FAA)  report  2  reviewing  its 
1978  airfield  pavement  design  guide  was  published.  This 
report  was  primarily  concerned  with  evaluating  the  existing 
FAA  pavement  design  procedures  for  high  traffic  volume 
pavements  (pavements  receiving  over  100,000  annual 
departures) .  In  the  design  review  and  data  collection 
phases  of  this  project,  it  was  pointed  out  that  there  was  no 
uniform  method  employed  in  the  design  of  the  pavement 
sections  with  regard  to  the  thickness  determination  of  the 
stabilized  base  layers.  Several  pavement  sections  were 
designed  using  detailed  theoretical  analysis,  and  others 
were  designed  using  empirical  methods  to  determine  the 
thickness  of  the  stabilized  base  system.  Test  sections  were 
also  used  in  some  designs  to  verify  the  assumptions  of  the 
theoretical  analysis. 

Most  of  these  facilities  were  constructed  in  the  early 
1970 's.  Since  that  time,  the  FAA  has  published  Advisory 
Circular  150/5320-6C,  dated  December  7,  1978.  This  circular 
includes  a  requirement  for  stabilized  bases  (when  gross 
aircraft  weight  exceeds  100,000  lb)  and  presents  a 
methodology  to  adjust  the  K-value  of  the  subgrade  for  the 
addition  of  stabilized  base.  This  method  is  easily  applied 
to  determine  the  thickness  of  the  concrete  slab  required. 
However,  determination  of  the  stabilized  base  thickness  is 
left  to  an  iteration  process.  A  suggestion  is  presented  in 
the  design  guide  that  the  thickness  of  the  base  course  be  as 
thick  as  the  required  slab.  This  statement  has  often  been 
interpreted  in  an  incorrect  fashion  by  some  designers.  The 
statement  has  been  taken  literally  where  a  concrete 
thickness  will  be  determined  for  a  given  K-value  and  then  a 
base  course  thickness  of  the  required  slab  designed  into  the 
project.  This  is  done  without  the  next  iteration  on  the 
base  course  to  determine  a  new  K-value. 

Other  problems  identified  in  the  current  approach  are 
as  follows: 

a.  There  is  no  adjustment  for  the  quality  of 
stabilized  materials  either  in  terms  of  modulus  of 
elasticity  or  modulus  of  rupture. 

b.  The  table  adjusts  the  K-value  to  a  maximum  value 
of  500  pci  beyond  which  there  is  no  adjustment. 
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g.  The  chart  usad  to  adjust  the  K-value  of  the 

subgrade  is  currently  truncated  at  a  stabilized 
base  thickness  of  12  in.  This  is  not  sufficient 
for  high-  volume  airfield  pavements. 

Purpose 

The  purpose  of  this  study  was  to  develop  a  design 
procedure  for  the  selection  of  the  thickness  of  a  stabilized 
base  course  material  for  rigid  airfield  pavements.  The 
design  method  was  developed  to  account  for  the  materials 
properties  of  the  stabilized  base  course,  and  to  be 
completely  compatible  with  the  elastic  layer  analysis  being 
used  by  the  US  Amy  Corps  of  Engineers. 

Soopo 

The  methodology  employed  in  this  study  is  applicable  to 
the  design  of  rigid  (portland  cement  concrete)  pavements 
utilizing  cement  treated  aggregate  as  the  base  course 
material.  The  program  developed  has  been  designed  using  a 
single-wheel  load,  a  dual-gear  load  (Boeing  727) ,  and  a 
dual-tandem  gear  load  (Boeing  767).  The  use  of  a  mixed 
traffic  equation  will  allow  the  procedure  to  be  used  in 
mixed  traffic  conditions  as  long  as  one  of  the  above  wheel 
loads  is  used  as  the  design  wheel  load. 

The  procedure  utilizes  elastic  layer  analysis  as  the 
theoretical  basis  of  the  method.  Thus,  it  is  feasible  that 
the  method  could  be  extended  to  other  stabilized  materials 
provided  that  their  materials  properties  are  within  the 
boundaries  used  in  this  study  and  that  their  primary  failure 
mechanism  is  tensile  fatigue. 

Organisation 

Section  2  presents  the  design  considerations  including 
the  overall  approach  to  the  problem,  selection  of  the  design 
factorial  and  selection  of  the  stress  analysis  procedure. 
Section  3  presents  the  results  of  the  stress  analysis  and 
discusses  the  relationships  of  the  variables.  The  stress 
determination  procedures  are  presented  in  Section  4. 

Sections  5  and  6  present  the  development  of  the  design 
procedure,  the  computer  programs,  and  the  comparison  of  the 
new  method  to  existing  design  approaches.  Section  7 
summarizes  the  findings  and  presents  conclusions  and 
recommendations  for  further  study. 
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DESIGN  CONSIDERATIONS 


Introduction 

In  order  to  develop  a  design  procedure  that  will  allow 
for  varying  the  stabilized  materials  properties,  the 
relevant  constitutive  properties  and  failure  mode  of  the 
materials  need  to  be  identified.  In  the  case  of  cement 
treated  bases,  Williams  1  has  identified  the  modulus  of 
elasticity,  modulus  of  rupture,  and  flexural  fatigue  as  the 
critical  materials  parameters  to  be  utilized  in  design. 
Other  researchers  have  identified  flexural  strain  (Mitchell 
and  Chen,.  Pretorios,  Raad)  ,0  as  the  failure  parameter. 
Costigan  1  :.n  his  work  for  the  US  Air  Force  used  the 
flexural  stress  approach.  In  order  to  maintain  consistency 
with  the  current  FAA  and  US  Army  Corps  of  Engineers  (CE) 
concrete  design  procedures,  the  flexural  stress  approach 
was  selected  for  use  in  this  study. 

The  basic  elements  of  this  study  are  as  follows: 

a.  Identify  variables  necessary  to  cover 
the  possible  range  of  materials 
properties. 

b.  Select  aircraft  type  to  be  used  in  the 
study . 

c.  Analyze  the  stress  state  for  the 
selected  aircraft  and  materials 
combinations. 

d.  Develop  the  design  procedure  based  on 
the  stress  state  of  the  paving  layers. 

e.  Compare  new  procedure  to  the  existing 
design  scheme. 

Thus,  the  major  parameters  that  must  be  considered  in 
the  design  process  are: 

a„  Thickness  of  the  concrete 

b.  Elastic  modulus  of  the  concrete 

c.  Thickness  of  the  cement  treated  base  (CTB) 

<L.  Elastic  modulus  of  the  CTB 

e.  Elastic  modulus  of  the  subgrade 
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f .  Traffic  parameters  (gear  type  and  load) . 

The  practical  range  of  each  of  these  factors  is 
required  to  develop  a  manageable  analysis  of  the?  tensile 
stress  conditions  for  combinations  of  concrete  surface  and 
CTB  base  course  layers. 

Factorial  Selection 

A  review  of  the  thicknesses  of  the  pavements  studied  in 
the  high  volume  traffic  design  project  was  performed  along 
with  the  associated  traffic  and  materials  properties  data. 

In  addition  to  this  review,  several  pavement  designs  were 
performed  to  determine  typical  concrete  thickness  to  define 
the  boundaries  of  the  problem.  This  synthesis  resulted  in 
the  following  parameters  to  be  included  in  the  factorial: 


£jl 

Thickness  of  the  concrete  (Hconc) ,  in. , 
22,  and  26. 

8,  10 

,  14, 

hj. 

Elastic  modulus  of  the  concrete  (Ec) ,  psi, 
4,000,000. 

SL i. 

Thickness  of  the  CTB  (Hctb) #  in.,  4,  8, 
and  22. 

10, 

14, 

sL. 

Elastic  modulus  of  the  CTB  (Ectb) t  ksi. 
1,000,  and  2,500. 

,  250, 

500, 

fijL 

Elastic  modulus  of  the  subgrade  (ESUb) , 
5,000,  10,000,  20,000,  and  30,000. 

psi. 

Gear  types,  single  wheel,  B-727,  B-767 

a*. 

Frictional  Factor  of  300  between  concrete  and 
and  CTB  and  subgrade. 

CTB 

This  resulted  in  a  total  of  1,200  cases  to  analyze  the 
stress  state  in  the  concrete  and  CTB  systems. 

Selection  of  Stress  Analysis  Procedure 

The  selection  of  the  theoretical  analysis  procedure  to 
define  the  stress  state  in  the  multilayer  system  was  the 
next  decision  point.  For  a  rigid  pavement  system,  there  are 
several  analysis  procedures  available  such  as  Westergaard, 
finite  element  methods,  or  elastic  layer  theory  to  analyze 
the  stress  state  in  the  pavement  structure.  Of  these 
approaches,  only  the  Westergaard  and  finite  element 
procedures  can  be  used  to  determine  the  pavement  behavior 
caused  by  the  edge  load  conditions.  However,  these  methods 
have  not  been  adopted  to  the  multilayer  case  with  different 
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interface  friction  conditions.  Finite  element  methods  have 
been  developed  which  can  handle  a  base  course  layer; 
however,  the  interface  condition  is  preset  at  full  bond  or 
no  bond. 

Analysis  of  the  non-destructive  testing  (NDT)  data 
obtained  on  the  high  traffic  volume  pavement  project  has 
shown  that  basin  matching  program  BISDEF  provides  better 
estimates  of  the  field  deflection  basins  if  the  friction 
value  between  the  layers  was  varied.  This  fact  in 
combination  with  the  other  factors  stated  above  lead  to  the 
use  of  the  "BISAR1"7  elastic  layer  procedure  for  use  in  the 
stress  analysis.  BISAR  is  structured  to  model  frictional 
interface  conditions  between  layers  on  a  continuous  scale. 

The  disadvantage  in  using  the  BISAR  system  is  its 
inability  to  handle  the  edge  (or  joint)  load  condition.  The 
inability  to  handle  the  edge  load  condition  has  been 
considered  a  serious  drawback  to  the  elastic  layer  analysis. 
Parker  et  al.  6  have  overcome  this  through  correlation  of 
edge  and  interior  load  conditions.  This  study  uses  a 
slightly  different  approach  which  accomplishes  the  same 
result.  This  will  be  discussed  below  and  in  Section  6. 

Stress  Analysis 

As  previously  mentioned  the  current  FAA  design 
procedures  are  based  on  the  edge  stress  loading  condition. 
This  is  the  load  condition  created  when  a  gear  or  wheel  is 
placed  on  a  free  edge  of  the  slab.  In  the  design  procedure 
a  value  of  75  percent  of  the  free  edge  load  is  used  to 
account  for  the  load  transfer  between  slabs.  In  order  to 
relate  the  FAA  edge  load  design  to  the  BISAR  elastic  layer 
interior  load  condition,  it  was  necessary  to  determine  the 
interior  stress  condition  for  the  case  of  a  concrete  slab  on 
subgrade.  For  this  set  of  systems,  the  frictional  factor 
was  set  to  the  full  slip  condition,  and  the  range  of 
thicknesses  was  changed  to  8,  10,  14,  18,  22,  and  26  in. 
over  the  same  subgrade  conditions.  This  resulted  in  an 
additional  72  cases  to  be  analyzed.  The  selection  of  the 
full  slip  case  was  made  to  model  more  closely  the  slab  on  a 
liquid  subgrade  condition.  This  condition  also  provides  a 
measure  of  safety,  since  the  full  slip  condition  results  in 
greater  stress  values. 

In  each  of  the  analysis,  the  maximum  tensile  stress  in 
the  concrete  and  the  CTB  was  located.  The  maximum  was  also 
located  for  the  slab-on-grade  cases.  For  the  multiple 
wheel  load  cases,  the  search  for  the  maximum  tensile  stress 
was  performed  in  several  locations.  Figure  1  presents  the 
general  layer  structure  and  interface  conditions  used  in  the 
analysis,  and  Figure  2  illustrates  the  search  locations  of 
the  multiple  wheel  assemblies.  As  shown  in  the  figure,  a 
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rigid  subgrade  layer  was  also  used.  This  layer  is  input  so 
that  surface  deflections  more  closely  match  the  field 
conditions . 

The  Friction  factor  of  300  represents  a  friction 
between  the  layers  of  about  30  percent.  The  value  is  a 
dimensionless  value.  The  BISAR  program  allows  varying  the 
friction  factor  from  0  to  1000.  The  300  value  was  based  on 
the  back  calculation  of  moduli  for  the  pavements  studied  by 
Kohn  . 
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Figure  1.  Pavement  structure  and  layer  interface  conditions 
used  in  the  stress  analysis 
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Figure  2.  Search  locations  for  the  maximum  tensile 
stress  for  the  multiple  wheel  assemblies 
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RESULTS  OF  STRESS  ANALYSIS 


The  results  of  the  1,272  systems  analyzed  are 
presented  and  discussed  in  this  section.  As  shown 
previously  in  Figure  2,  multiple  search  locations  were  used 
in  the  analysis.  The  pavement  configuration  and  interface 
condition  were  presented  in  Figure  1. 

The  resulting  maximum  tensile  stresses  are  presented 
in  Appendix  A.  Table  A1  is  illustrative  of  the  set  of 
results.  This  table,  presented  in  the  appendix,  includes 
the  results  for  the  pavements  containing  the  CTB  layer. 

All  stresses  presented  are  the  maximum  radial  tensile 
stress  for  the  condition  represented. 

In  order  to  initiate  the  development  of  the  design 
process,  it  was  already  determined  that  the  maximum  tensile 
stress  of  each  layer  would  be  the  major  design  factor. 

Prior  to  starting  the  analysis  procedures,  graphic 
representations  of  the  stress  in  the  concrete  (STRC)  and 
stress  in  the  base  course  (STRB)  were  prepared.  Plots  were 
developed  of  both  the  stress  factors  against  each  of  the 
major  variables.  A  complete  set  of  graphs  is  presented  in 
Appendix  B.  Figures  3  to  22  represent  the  stress  condition 
for  the  single  wheel  load.  Figures  3  to  18  represent  the 
mean  pavement  parameters  used  in  the  stress  factorial,  and 
Figures  19  to  22  present  the  conditions  at  the  boundaries 
of  the  factorial. 

Stress  in  the  Concrete 

Figures  3  to  10  present  the  primary  relationships  for 
the  stress  in  the  concrete  slab.  As  the  curves  in  the 
figures  show,  the  relationship  between  STRC  and  Hconc  is 
certainly  a  power  function.  From  the  theoretical  analysis 
of  the  condition  it  is  a  function  of  VHZconc.  This  can  be 
determined  through  the  review  of  the  Westergaard  equations. 
Scott  presents  an  excellent  summary  of  these  equations 
and  other  analytical  approaches  to  calculating  the  tensile 
stress  in  a  slab  resting  on  either  a  Winkler  (liquid)  or 
elastic  foundation. 

Some  other  general  observations  can  be  made  regarding 
the  stress  in  the  concrete  when  viewing  the  plots.  The 
elastic  modulus  of  the  subgrade  (E^)  alone  does  not  have  a 
major  effect  on  the  stress  (Figure  9) .  However,  there  is 
an  interaction  effect  of  Esujj  and  the  elastic  modulus  of  the 
base  course  (E^) .  The  thickness  of  the  stabilized  base 
has  a  fairly  linear  reduction  effect  on  the  value  STRC,  and 
the  rate  of  the  reduction  is  dependent  on  the  value 
(Figure  4) . 
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Figure  5.  Relationship  between  CTB  moduli  and  stress  in 
the  concrete  slab  for  all  concrete  thicknesses 
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Figure  6.  Relationship  between  CTB  moduli  and  stress  in 
the  concrete  slab  for  all  CTB  thicknesses 
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Figure  15.  Relationship  between  CTB  thickness  and  stress 
in  the  CTB  for  all  subgrade  moduli 
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Figure  16.  Relationship  between  subgrade  moduli  and 
stress  in  the  CTB  for  all  CTB  moduli 
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Figure  19.  Relationship  between  concrete  thickness  and  stress 
in  the  concrete  slab  for  all  subgrade  moduli 
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evaluated  at  the  boundaries  of  the  factorial 
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Figure  22. 
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All  these  relationships  are  generally  diminished  as  the 
slab  thickness  becomes  very  large  (Figure  5) .  Thus,  for  a 
given  load  condition,  the  concrete  slab  thickness  can 
overpower  all  parameters  if  it  is  extremely  thick.  This 
would  happen  in  the  case  of  over-designed  pavements. 

An  interesting  result  is  found  when  reviewing  the 
graphs  for  the  dual  and  dual  tandem  gears.  It  can  be  seen 
that  the  same  general  trends  are  produced  for  the  multiple- 
wheel  loads.  This  was  somewhat  expected  since  it  was  shown 
by  Woinowsky-Krieger  12  that  the  tensile  stress  directly 
under  the  load  plate  could  be  estimated  by  an  equation  only 
involving  the  thickness  of  the  slab  and  no  other  parameters. 
However,  the  effect  of  overlapping  stresses  from  the 
multiple  wheels  was  unknown  prior  to  the  analysis. 

Stress  in  the  Base  Course 

Figures  11  to  18  present  the  same  primary  relationships 
as  Figures  3  to  10  for  the  stress  in  the  STRB .  In  general, 
similar  relationships  hold.  The  thickness  of  the  overlying 
slab  has  a  large  influence  on  the  value  of  the  STRB.  The 
stress  in  the  base  as  a  function  of  the  thickness  of  the 
base  was  found  to  be  the  most  difficult  relationship.  As 
shown  in  Figures  12  and  17,  this  relationship  tends  toward  a 
parabolic  form.  Thus,  transformation  to  a  linear  function 
would  be  difficult.  However,  in  the  range  of  the  factorial, 
this  nonlinear  behavior  is  not  fully  developed. 

Unlike  the  value  STRC,  the  stress  in  the  base  is  not 
influenced  greatly  by  the  thickness  of  the  base  until  better 
quality  (higher  Ectfc>)  materials  are  used.  Figures  12  and  14 
illustrate  this  behavior. 

Again,  the  relationships  for  the  dual  and  dual  tandem 
gears  were  found  to  follow  the  same  general  trends.  The 
slopes  and  magnitudes  of  the  changes  in  stress  are  different 
but  the  shapes  do  remain  fairly  constant. 

Boundary  Relationships 

The  previous  set  of  figures  were  developed  for  the  mean 
pavement  parameters.  However,  before  pursuing  any  analysis 
procedures,  conditions  representing  the  boundaries  of  the 
factorial  were  developed.  These  relationships  are  presented 
in  Figures  19  to  22.  These  graphs  represent  the 
relationships  between  the  major  variables  and  the  stresses 
for  the  conditions  at  the  four  corners  of  the  factorial.  As 
shown,  the  relationships  hold  in  a  general  sense,  but  it  is 
shown  (Figure  20)  that  interaction  effects  are  present. 

That  is  relationships  which  were  once  approximately  parallel 
over  the  range  of  thicknesses  begin  to  cross  and  become  non- 
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parallel.  Also,  the  nonlinear  behavior  of  the  STRB  value 
becomes  magnified  at  the  edges  of  the  factorial. 

Summary 

These  graphics  present  the  general  relationship  between 
the  stress  values  (STRC  and  STRB)  and  the  materials  and  slab 
characteristics  governing  the  design  problem.  The  review  of 
these  graphs  lead  to  the  possibility  that  a  single  form 
could  be  used  for  determining  the  stress  values  for  any  gear 
configuration  (the  equation  was  modified  by  a  gear  factor) . 
This  can  be  expressed  in  the  following  form: 

Stress  =  f (G) *f (h, Esub* Ectb) 


where  f(G)  =  gear  factor 
f (h/EgubfEctb  )  =  general  function  in  terms 

of  materials  properties. 

The  other  conclusion  was  that  the  determination  of  the 
stress  in  the  base  would  be  a  more  complex  function  since 
the  parameter  exhibited  the  parabolic  nature  especially  at 
the  extremities  of  the  factorial. 

At  this  point  in  the  project  it  was  determined  that  the 
two  primary  methods  to  consider  were  an  analytical  or  a 
statistical  approach.  Graphic  procedures  were  put  aside  at 
this  point  due  to  the  complex  interactive  effects  of  the 
variables. 
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STRESS  DETERMINATION  PROCEDURES 


Introduction 

The  objective  is  to  design  a  two-layer  rigid  (surface 
and  base)  pavement  system  that  will  provide  equivalent 
performance  to  a  one-rigid  layer  (surface)  pavement  system 
as  determined  by  the  current  FAA  procedures.  The  new  design 
method  is  based  on  obtaining  equal  maximum  tensile  stress  at 
the  bottom  of  the  concrete  layer  in  a  two-layer  pavement 
system  without  exceeding  the  allowable  tensile  stress  levels 
in  the  rigid  base. 

This  section  describes  the  methods  investigated  and  the 
method  selected  to  predict  the  maximum  tensile  stress  in 
each  of  the  pavement  layers. 

Methodology 

Initially,  two  theoretical  approaches  were  investigated 
in  order  to  predict  the  required  pavement  behavior.  The  two 
approaches  were  both  an  equivalent  thickness  approach.  One 
method  relied  on  an  equal  rigidity  concept  and  the  other  was 
based  on  an  equivalent  stress  concept.  These  methods  are 
used  frequently  in  overlay  design  to  determine  the  thickness 
of  the  overlay  slab.  In  many  respects  the  design  of  a  two 
rigid  layer  system  is  very  similar  to  overlay  design.  It 
was  thought  that  modification  of  the  equivalent  thickness 
for  gear  load  and  subgrade  effects  could  be  developed  to 
provide  a  system  of  equal  stress  and  therefore  equal 
performance.  The  equation  used  in  this  analysis  is  as 
follows: 

He  =  (Hsn+Hctbn)  *  F(Eb)  (1) 

where 

He  =  the  equivalent  thickness 

Hs  =  the  surface  thickness 
Hctb  =  the  base  thickness 

n  =  a  power  determined  by  the  method  used 

F(Eb)  =  a  modular  ratio  term. 

This  equation  can  be  developed  from  simple  beam  theory 
as  presented  by  Kohn  and  Rollings  .  Once  the  equivalent 
thickness  was  determined,  the  equivalent  pavement  system 
were  treated  as  a  linear  elastic  slab  supported  by  a  dense 
liquid.  Westergaard 1 s  stress  equation  is  then  used  to 
determined  the  required  pavement  behavior. 
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Neither  approach  was  successful  in  predicting  the 
maximum  tensile  stress  at  the  bottom  of  the  concrete  slab  as 
determined  by  the  elastic  theory.  The  major  problem  in 
developing  the  necessary  modifying  functions  was  the 
interaction  effects  at  the  boundary  of  factorial.  Thus,  the 
linear  regression  analysis  was  selected  to  predict  the 
tensile  stress  with  an  acceptable  level  of  accuracy. 

Regression  Analysis 

Five  factors  were  evaluated  in  this  study  to 
investigate  their  influence  on  the  maximum  tensile  stress  at 
the  bottom  of  the  concrete  layer.  The  maximum  tensile 
stress  has  been  clearly  identified  as  having  a  great 
influence  cn  the  structural  performance  of  rigid  pavement 
systems.  The  five  variables  are: 

a.  Loading  or  gear  type,  GT 

b.  Modulus  of  elasticity  of  the  CTB,  Ectb 

c.  Thickness  of  the  concrete  surface,  Hconc 

d.  Thickness  of  the  CTB,  Hctb 

e.  Elastic  modulus  of  the  subgrade,  ESU£, 

The  linear  regression  analysis  was  performed  using  the 
Statistical  Package  for  Social  Studies  (SPSS)  microcomputer 
program.  SPSS  is  a  widely  accepted  program  used  in  both  the 
social  science  and  engineering  fields  for  statistical 
analysis.  The  program  offers  a  wide  range  of  regression 
methods,  and  allows  for  building  transformation  variables. 
The  steps  involved  in  the  selection  of  the  final  form  of  the 
equation  were  as  follows: 

a.  Identification  of  a  large  set  of  variables  which 
are  thought  to  influence  the  stress  state  of  the 
pavement  system. 

b.  Selection  of  the  variables  which  are  thought  to 
have  the  most  influence  on  the  pavement  response. 
Three  factors  controlled  the  selection  variables: 

(1) .  variables  should  be  a  function  of  pavement 

geometry  and/or  material  properties. 

(2)  .  the  form  of  the  equation  should  be  the  same 

for  all  gear  types. 


(3) .  the  resulting  equation  should  provide  the 
least  variability  from  the  stresses 
generated  in  the  factorial  design  for  each 
gear  type. 

c.  Once  the  equation  is  developed,  a  statistical 
model  and  adequacy  check  is  performed.  If  the 
model  is  statistically  inadequate,  return  to 
step  b. 


Selection  of  Variables 


Based  on  the  figures  shown  in  the  previous  section,  it 
was  obvious  that  strong  interactions  exist  between  the 
factors  considered  in  this  study.  The  need  to  consider 
interactions  between  the  factors  was  also  identified  by 
obvious  patterns  in  early  stage  residual  plots  and  by  strong 
intercorrelation  between  the  factors  in  the  correlation 
matrix  when  only  the  main  effects  were  considered.  The 
figures  also  illustrated  that  power  and  log  functions 
represent  most  relationships  between  the  different  factors. 

Based  on  these  observations,  a  large  set  (over  100 
initially)  of  variables  consisting  of  interactions  between 
the  factors  investigated  and  other  related  variables 
created.  Trial  and  error  procedures  were  used  to  select  the 
appropriate  power  functions.  Factors  that  were  observed  in 
the  equivalent  thickness  approach  were  also  used  to  develop 
variables. 

Westergaard ' s  stress  solution  for  an  interior  load  on  a 
slab  supported  by  a  dense  liquid  indicated  that  the  radius 
of  relative  stiffness  of  the  slab  is  a  very  important 
parameter.  The  radius  of  relative  stiffness  3  the  inverse 
of  the  parameter  lambda  determined  from  the  displacement 
equation  of  a  beam  on  a  Winkler  foundation.  The  unit  of  the 
radius  of  relative  stiffness  is  length.  It  is  a  measure  of 
the  interaction  effect  between  the  slab  and  the  subgrade. 

If  the  distance  is  large  the  deflection  of  the  soil  and  beam 
will  be  over  a  large  area,  whereas  if  the  distance  is  small 
the  influence  of  the  load  will  be  localized.  In  the  case  of 
the  three  layer  system  the  value  of  the  radius  of  relative 
stiffness  calculated  using  the  properties  of  the  surface  and 
the  subgrade  is  not  sufficient  to  determine  the  stress  in 
both  layers. 

An  attempt  was  made  to  develop  an  equivalent  radius  of 
relative  stiffness  for  the  two  rigid  layers  sing  the 
equivalent  thicknesses  described  earlier.  However,  this 
attempt  was  unsuccessful  since  it  did  not  account  for  the 


variability  observed  in  the  pavement  response  (STRC) .  A 
radius  of  relative  stiffness  was  determined  for  the  surface 
(L})  and  base  (L2)  slabs  assuming  they  are  supported  by  the 
same  dense  liquid.  Interactions  between  Lj  and  L2  or 
between  functions  of  Li,  L2  and/or  other  variables 
accounted  for  most  variability  in  STRC. 

The  equivalent  slab  thickness,  based  on  equal  tensile 
stress  at  the  bottom  of  the  slab,  indicated  that  the  ratio 
of  the  modulus  of  the  different  layers  is  an  important 
variable.  Consequently,  these  ratios  were  considered  in 
the  analysis. 

The  second  step  in  the  variable  selection  process 
involved  using  the  stepwise  variable  elimination  procedure 
provided  by  SPSS.  This  procedure  reduced  the  number  of 
variables  to  a  manageable  number  without  greatly  affecting 
the  accuracy  in  predicting  STRC. 

The  third  step  was  to  check  the  adequacy  of  the  model 
to  ensure  the  selection  of  a  sound  statistical  model.  This 
part  of  the  process  is  required  to  eliminate  gross  errors 
in  the  prediction  of  the  stress.  Normal  probability  plots 
and  residual  plots  were  used  to  complete  this  part  of  the 
process. 

The  main  components  of  the  variables  that  entered  the 
equations  are  listed  below: 


Er^  =  ratio  of  concrete  modulus  to  the  CTB 
modulus 


Er2 


Hc  an<*  Hctb 


L 


L1'L2 


K 


ratio  of  CTB  modulus  to  elastic  modulus 
of  the  subgrade 

thickness  of  concrete  slab  and 
stabilized  base,  respectively 

radius  of  relative  stiffness  of  the 
concrete  slab  for  the  one  rigid  layer 
pavement  system 

radius  of  relative  stiffness  of  the 
concrete  slab  and  stabilized  base, 
respectively,  for  the  two  rigid  layer 
pavements  systems 

modulus  of  subgrade  reaction 


Esub  =  Elast-i-c  modulus  of  the  subgrade 


Results  of  the  One-Rigid  Layer  System 


The  maximum  tensile  stress  in  the  concrete  slab  was 
predicted  for  this  case.  One  equation  form  was  developed 
for  all  three  gear  configurations.  Naturally,  different 
regression  coefficients  were  computed  for  each  gear  type. 
The  general  form  of  the  equation  is  as  follows: 

SC  «  P/Pl  {CONSTANT  +  al*ln(Esub)  +  l/Hconc2* [a3*L  +  a4*K] 

+  ln(L)*[a5*Hconc+  a2*ln(Esub)]  +  ae.l/H^3)  (2) 

Where 

L  -  (  tEconc  *  Hconc3l  /  U*  d  '  “2)  *  *1  )025 
log  K  -  (log  Esub  -  1.415)  /  1.284 

ai  -  regression  coefficients 

SC  —  maximum  tensile  stress  slab-on-grade 

Er  =*  Modulus  ratio 

Hconc  *  concrete  thickness 

L  =  radius  of  relative  stiffness 

K  =  modulus  of  subgrade  reaction 

P  =  design  gross  aircraft  load 

PI  =  84,000  lb  for  single  wheel  aircraft 
190,000  lb  for  dual  wheel  aircraft 
305,000  lb  for  dual  tandem  aircraft 


The  values  of  the  regression  coefficients  ai,  the 
correlation  coefficients  R  ,  and  the  standard  error  for  each 
gear  type  are  given  in  Table  1.  The  regression  equations 
provided  excellent  predictions  for  the  maximum  tensile  stress 
as  seen  ir.  the  values  of  RZ  and  standard  error.  A  standard 
error  of  2.7  psi  was  computed  which  represents  0.63  percent 
Of  the  mean  stress  value  for  the  dual-wheel  gear.  The 
maximum  standard  error  calculated  was  equal  to  0.75  percent 
and  corresponded  to  the  single  wheel  gear  type.  Plots 
illustrating  the  agreement  between  predicted  and  computed 
stress  values  are  shown  in  Figure  23  through  25  for 
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the  single  wheel,  dual  wheel,  and  tandem  gear  types, 
respectively.  Additional  plots  illustrating  variance  in  the 
residuals  of  the  selected  models  are  shown  in  Figures  26 
through  28. 


Table  1 


Regression  Coefficients.  Interior  Stress  for  Slab-On-Grade 


GEAR  TYPE 

REGRESSION 

SINGLE 

DUAL 

DUAL 

COEFFICIENTS 

WHEEL 

WHEELS 

TANDEM 

al 

79.87059 

70.87699 

-104.54994 

a2 

-20.83124 

-13.20461 

41.08126 

a3 

1043.40315 

1886.23113 

1754.38235 

a4 

-28.99054 

-35.41080 

31 . 44690 

a5 

1.11244 

0.91015 

-1.86926 

a6 

33382.69811 

-74244.45618 

Constant 

-54.80360 

-307.15867 

-601.92659 

R2 

0.99988 

0.99991 

0.99987 

Standard 

Error, 

psi 

2.08718 

2.70057 

2.09560 

Standard 

Error, 

%  error  0.75 

0.63 

0.67 

It  is  important  to  note  that  no  strong  patterns  were 
observed  in  the  residual  plots.  This  indicates  that  the 
selected  model  is  adequate  to  predict  the  tensile  stress  in 
the  concrete  slab  when  the  design  inputs  are  within  the 
valid  range  of  the  investigated  factors.  As  indicated 
earlier  in  this  report,  the  range  of  design  inputs  was 
selected  to  include  most  practical  pavement  design 
applications . 
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REGRESSION  EQUATION 

Figure  23.  Comparison  of  predicted  versus  actual  stress  values  for  the 
case  of  a  single-wheel  gear  and  the  slab-on-grade  condition 
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REGRESSION  EQUATION 

Figure  24.  Comparison  of  predicted  versus  actual  stress  values  for  the 
case  of  a  dual-wheel  gear  and  the  slab-on-grade  condition 
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REGRESSION  EQUATION 

Figure  26.  Residual  plot  obtained  from  the  statistical  model  selected 
for  the  single-wheel  gear  (slab-on-grade  condition) 


Results  of  the  Two-Rigid  Layer  Pavement  system 


For  this  case/  the  maximum  tensile  stress  in  the 
concrete  slab  and  the  cement-treated  base  were  predicted. 
The  equation  was  selected  to  provide  the  best  fit  using  one 
equation  form  for  all  gear  types. 

The  equation  for  the  STRC  is  as  follows: 

STRC  *  P/P CONSTANT  +  al  InfL^/Lj  +  H^^lnfL})  {  a2  +  a5 

2 

*  L2]  +  l/Hconc  *  {  a3*  Lj  +  a6*  Er j*Lj  +  a9*  L2  + 
a^i*  K*(Lj/L2)  t  a8*  L2*L j  }  +  B  *  {  a4  +  a7*  Er2+ 

a10*  l2  +  a12*  (Erl/Esub>°'2S  >>  <3> 


where 

ai  =  regression  coefficients  (see  Table  2) 

Erl  =  Econc  /  Ectb  ;  Er2=  Ectb  /  Esub 
L1  -  <  tEconc  *  "con/l  '  I12  <2  ‘  “2>  *  K1  >°'25 

l2  “  (  tEctb  *  Hctb3>  /  12  <l  *  u2>*kJ  >a25 

log  K  =  (log  Esub  -  1.415)  /  1.284] 

B  -  ^"conc0'33  *  ln  (VW  *  <HCtl/Hconc)  2*  (Er ,)  "25 
STRC  =  maximum  tensile  stress  in  the  concrete 


HCOnc  =  concrete  thickness 

Hctb  =  cement  treated 

L  =  radius  of  relative  stiffness 


K  =  modulus  of  subgrade  reaction 

P  =  design  gross  aircraft  load 

Pj  =  84,000  lb  for  single  wheel  aircraft 

190,000  lb  for  dual  wheel  aircraft 
305,000  lb  for  dual  tandem  aircraft 
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The  values  of  the  regression  coefficients  along  with 
other  statistical  parameters  are  presented  in  Table  2. 
Relatively  good  correlation  was  obtained  for  all  gear  types 
considering  the  complexity  of  the  three  layer  pavement 
system.  Minor  changes  in  any  of  the  parameters 
representing  the  different  pavement  layers  causes 
significant  changes  in  the  observed  pavement  response.  The 
maximum  standard  error  computed  for  the  regression 
equations  was  equal  to  2.55  percent  of  the  mean  tensile 
stress.  This  value  corresponds  to  the  dual-tandem  gear. 

Table  2 

Regression  Coefficients.  Interior  Stress  in  the 
Concrete  Slab  for  the  Two-Rigid  Laver  System 


REGRESSION 

COEFFICIENTS 

GEAR  TYPE 

SINGLE 

WHEEL 

DUAL 

WHEELS 

DUAL 

TANDEM 

al 

3152.60047 

3903.81004 

5463.62830 

a2 

0.74036 

1.06777 

a3 

1206.42011 

1839.96268 

1304.04836 

a4 

40.117044 

79.37649 

66.63757 

a5 

2 . 891281E 

-03 

a6 

3.8000969811 

6.99535 

5.63003 

a7 

0.02337 

0.03171 

0.03013 

a8 

-5.30015 

a9 

-925.25867 

-1421.66496 

-820.59149 

alO 

0.45975 

0.73535 

0.36554 

all 

-5.58656 

-10.42594 

-7.26875 

al2 

-44.92398 

-195.49499 

-139.49406 

Constant 

-307.75101 

-388.66280 

-41.35524 

R2 

0.99912 

0.99905 

0.99814 

Standard  Error,  psi 

3.66431 

5.55339 

5.26719 

Standard  Error,  %  error  1.96 

1.98 

2.55 

The  high  degree  of  variability  associated  with  this 
gear  type  can  be  justified  by  considering  the  following 
factors.  The  predicted  pavement  response  is  the  maximum 
tensile  stress  in  the  concrete  slab  and  not  the  tensile 
stress  at  one  location  within  the  pavement  system.  As  the 
thicknesses  and  material  properties  changed  in  the  factorial 
design,  the  location  of  the  maximum  tensile  stress  within 
each  of  the  rigid  layers  also  changed.  It  is  believed  that 
this  variation  in  the  maximum  tensile  stress  location 
accounts  for  most  of  the  variability  observed  in  the 
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regression  results.  The  same  phenomenon  was  also  observed, 
but  to  a  lesser  degree  in  the  dual  wheel  results.  The 
magnitude  of  the  variance  caused  by  this  phenomenon  is  much 
greater  for  the  stress  in  the  base. 

The  largest  magnitude  of  the  standard  error  observed 
for  all  gear  types  is  equal  to  5.55  psi  which  is  acceptable 
for  all  practical  design  applications. 

Figures  29  through  31  show  the  agreement  between  the 
predicted  and  computed  stress  values  for  all  three  gear 
types.  Figures  32  through  34  present  the  residuals  for  each 
of  the  three  equations.  These  residuals  do  not  show  any 
gross  violations  for  the  assumptions  (i.e.  normality  and 
equal  variance  of  the  error  term)  made  during  the 
statistical  analysis. 

The  stress  in  the  cement-treated  base  course  was  also 
predicted.  For  this  value  the  natural  log  of  the  maximum 
tensile  stress  in  the  base  was  predicted.  The  logarithmic 
transformation  was  required  to  reduce  the  variance  in  the 
error  term  which,  otherwise,  would  have  violated  one  of  the 
basic  assumptions  of  statistical  analysis.  This  high 
magnitude  of  the  variance  can  be  explained  in  part  by  the 
change  in  the  maximum  tensile  stress  which  was  addressed  in 
the  previous  paragraphs  and  the  nonlinear  behavior  of  the 
function  at  the  extremities  of  the  factorial. 

The  selected  model  fits  all  gear  types  with  an 
acceptable  level  of  accuracy.  The  regression  equation 
developed  is  as  follows: 

Ln(STRB)  -  CONSTANT  +  (Erx/Egub) 0 * 25  {a5  Hconc*ln(L1)  +  a6  * 
Er2  +  alO  *  Erj.  +  al2>  +  K/Hconc2  (al  Hctb2+  a4 
L!/L2}  +  Ln(Li)*  Hconc{a3  Li/Hconc2  +  a7  B  +  all 
LnCLiJ/Li)  +  Erx  *  {a2  +  a8  L2)+a9  Ln^)/!,!  (4) 


STRB=  P/P1(exp(Ln(STRB) )  (5) 


where 

Erl  =  E  Conc  /  E  ctb 
Er2  =  E  ctb  /  E  sub 
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L1  -  I  tEconc  *  JW3]  /  [12  (1  -  u2).  k]  l0'25 
1*2  -  (  [E^b  *  H^b3)  /  [12  <l-u2)*k]  )0-25 

log  k  -  (log  Esub  -  1.415)  /  1.284 

B  “  VH«c°'33  *  ln<VHconC>*<Hctl/Hconc)2*<Erl> 
ai  =  regression  coefficients  (see  Table  3) 

STRB  *  maximum  tensile  stress  in  the  CTB 

Hconc  *  concrete  thickness 

Hctb  =  cement-treated  base  thickness 

L  =  radius  of  relative  stiffness 

k  =  modulus  of  subgrade  reaction 

P  *  design  gross  aircraft  load 

Pj>=  84,000  lb  for  single-wheel  aircraft 
190,000  lb  for  dual-wheel  aircraft 
305,000  lb  for  dual-tandem  aircraft 


-.25 
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EAST  SQUARE  ESTIMATES  FOR 
DUAL  GEAR 


RESIDUALS  OF  LEAST  SQUARE  ESTIMATE 
FOR  SINGLE  WHEEL 


sivnais3d 
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REGRESSION  EQUATION 

Figure  32.  Residual  plot  obtained  from  the  statistical  model  selected 
for  the  single-wheel  gear  (two-rigid  layer  system) 


RESIDUALS  OF  LEAST  SQUARE  ESTIMATE 
FOR  DUAL  GEAR 


sivnais3a 
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REGRESSION  EQUATION 

Figure  33.  Residual  plot  obtained  from  the  statistical  model  selected 
for  the  dual-wheel  gear  (two-rigid  layer  system) 


RESIDUALS  OF  LEAST  SQUARE  ESTIMATE 
FOR  DUAL  TANDEM  GEAR 


sivnais3u 
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REGRESSION  EQUATION 

Figure  34,  Residual  plot  obtained  from  the  statistical  model  selected 
for  the  dual-tandem  gear  (two-rigid  layer  system) 


Table  3 


Regression  Coefficients.  Interior  Stress  in  the  CTB 
for  the  Two-Riaid  Laver  System 


GEAR  TYPE 

REGRESSION 

SINGLE 

DUAL 

DUAL 

COEFFICIENTS 

WHEEL 

WHEELS 

TANDEM 

al 

-1.54962E-04 

-1. 14871E-04 

-9 . 69983E-05 

a2 

-.09001 

-0.08893 

-0.09306 

a3 

0.32334 

0.33870 

0.36498 

a4 

-0.02994 

-0.02350 

-0.02921 

a5 

-0.01850 

-0.01829 

-9 . 29258E-03 

a6 

-6 . 82371E-03 

-7 . 46144E-03 

-3 . 35922E-03 

a7 

-7 . 93621E-03 

-7 . 78525E-03 

-7 . 47100E-03 

a8 

2 . 778992E-04 

5 . 404066E-04 

7 . 70275E-04 

a9 

22.59014 

20.06895 

17.86815 

alO 

0.47821 

0.46020 

0.45719 

all 

-0.15486 

-0.12404 

-0.08699 

al2 

-22.51655 

1.92671 

-23.879356 

Constant 

1.57794 

-388.66280 

1.46192 

R2 

0.99827 

0.99852 

0.99772 

Standard 

Error, 

0.03965 

0.03443 

0.03970 

psi 

Standard 

Error, 

%  error 

1.11 

0.84 

1.0 

The  maximum  standard  error  was  equal  to  1.11  percent  of 
the  logarithmic  mean  of  the  stress  values.  A  high  degree  of 
correlation  was  observed  between  the  predicted  and  computed 
stress  values  as  indicated  by  the  high  Rz  values.  Figures 
35  through  37  show  the  close  agreement  between  the  predicted 
and  computed  stresses.  Residual  plots  presented  in  Figures 
38  through  40  do  not  show  any  trends  but  indicate  a  small 
error  over  the  range  of  the  regression. 
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RESIDUALS  OF  LEAST  SQUARE  ESTIMATE 
FOR  SINGLE  WHEEL 
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RESIDUALS  OF  LEAST  SQUARE 
FOR  DUAL  GEAR 
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Figure  39.  Residual  plot  obtained  from  the  st< 
for  the  dual -wheel  gear  (two-rigid 


RESIDUALS  OF  LEAST  SQUARE  ESTIMATE 
FOR  DUAL  TANDEM  GEAR 


Figure  40.  Residual  plot  obtained  from  the  statistical  model  selected 
for  the  dual-tandem  gear  (two-rigid  layer  system) 


Load  Effacts 

The  factorial  design  did  not  include  the  aircraft  load 
as  one  of  the  major  factors.  Since  BISAR  was  the  selected 
behavior  model;  linear  elasticity  is  assumed.  However, 
problems  could  arise  through  the  interaction  of  two  factors. 
Partial  interface  slip  conditions  were  used  in  the  analysis. 
This  fact  alone  would  not  cause  nonlinear  behavior.  Coupled 
with  the  fact  that  the  maximum  tensile  stress  was  used 
could,  however,  develop  some  problems  with  regard  to 
positional  change  of  the  maximum  load.  If  these  factors  did 
cause  positional  shifts,  a  modified  equation  for  loading 
conditions  would  be  necessary.  To  validate  the  linearity 
assumption,  two  factorials  were  analyzed  using  BISAR — one 
for  the  concrete  (surface  layer)  and  the  other  for  the 
cement-treated  base.  Primarily,  the  boundaries  of  the 
original  factorial  were  searched.  These  areas  were  found  to 
contain  the  elements  of  positional  shifts  in  the  original 
computer  runs. 

The  results  of  these  runs  are  presented  in  Tables  4 
and  5.  The  ta/»  es  present  the  ratio  of  the  maximum  tensile 
stress  determined  at  the  full-  and  half-load  limits.  As 
shown,  the  ratio  of  the  stresses  calculated  for  the  full- 
and  half-load  case  is  nearly  constant  at  a  value  of  0.5. 
Thus,  the  assumption  of  linearity  is  not  interfered  with  due 
to  the  positional  and  friction  factor  considerations. 

Summary 

Although  the  selected  regression  equations  contain  a 
large  number  of  variables  to  reduce  the  error  in  prediction, 
most  of  the  variables  are  used  in  current  design  procedures 
and  are  easy  to  compute.  In  addition,  a  microcomputer 
program  will  be  provided  to  solve  efficiently  the  regression 
equations  and  to  select  the  required  base/surface 
thicknesses  for  a  particular  foundation  and  loading 
condition. 

The  errors  associated  with  the  regression  equations  are 
judged  to  be  acceptable  for  all  practical  design  purposes. 
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ess  Comparison  at  Load  =  P  and  Load  =  p/2  (Conor 
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Table  5 


Stress  Comparison  at  Load  «  P  and  Load  =  P/2  ( CTB1 


E  SUb  .  ksl 

S 

30 

ECTB.  hsi 

260 

2600 

HC,  in. 

8 

26 

8 

26 

GEAR 

Typg 

HOTS,  In. 

LOAD 

IVP£ 

4 

22 

4 

22 

P/2 

46.900 

6.600 

167.000 

22.700 

DUAL 

P 

93.700 

13.000 

314.000 

45.400 

STRESS  RATIO 

0.501 

0.500 

0.500 

0.500 

P/2 

34.800 

8.000 

1  OS.  000 

22.200 

DUAL 

TANDUM 

P 

69.600 

16.100 

217.000 

44.500 

STRESS  RATIO 

0.600 

0.499 

0.490 

0  499 

DEVELOPMENT  07  THE  DESIGN  PROCEDURE 


This  section  describes  the  assumptions  made,  the 
limitations,  and  the  steps  involved  in  the  new  design 
procedure.  This  section  will  also  include  descriptions  of 
two  microcomputer  programs  provided  to  perform  the  design  of 
new  rigid  pavements  and/or  the  evaluation  of  existing 
pavements . 

The  new  design  method  was  developed  to  supplement 
existing  rigid  pavement  design  methods  used  by  the  F\A  and 
not  to  replace  them.  The  new  method  is  based  on  providing  a 
two  rigid  layer  pavement  structure  equivalent  in  terms  of 
performance  to  the  concrete  slab  (one  rigid  layer)  designed 
using  current  FAA  procedures.  The  equivalence  between  the 
two  pavement  structures  is  achieved  by  equating  the  maximum 
tensile  stress  in  the  concrete  slab  without  exceeding  the 
tensile  strength  of  the  CTB. 

Assumptions 

Current  rigid  pavement  design  procedures  are  based  on 
limiting  the  maximum  tensile  stress  at  the  edge  of  the 
pavement.  The  critical  stress  condition  occurs  when  the 
load  is  located  as  close  as  possible  to  the  edge  of  the 
pavement.  However,  because  BISAR  uses  elastic  theory  to 
compute  the  pavement's  behaviour  to  loading,  the  edge 
loading  condition  could  not  be  modeled.  This  limitation 
arises  from  an  assumption  made  by  the  elastic  theory, 
namely,  the  pavement  dimensions  extend  to  infinity  in  the 
horizontal  plane.  For  this  reason,  interior  stresses  were 
evaluated  instead  of  edge  stresses.  The  main  assumption 
made  during  the  course  of  developing  this  design  procedure 
is  that  similar  performance  will  be  achieved  by  equating  the 
maximum  interior  tensile  stress  at  the  bottom  of  the 
concrete  layer  in  both  pavement  structures.  If  the  maximum 
interior  tensile  stress  is  equal  in  both  pavement 
structures,  the  maximum  edge  tensile  stress  is  also  equal  in 
both  pavement  structures. 

Another  limitation  to  this  design  procedure  is  that  the 
method  is  valid  only  within  the  range  of  the  factorial 
design.  This  limitation  arises  from  the  use  of  regression 
analysis  to  predict  the  required  pavement  response. 

However,  the  same  approach  can  be  used  if  the  stress  can  be 
accurately  predicted  outside  the  range  of  the  factorial 
design  using  such  methods  as  elastic  theory. 

Strength  and  failure  criteria  of  the  CTB  have  been 
obtained  from  limited  published  data.  These  relationships 
should  be  modified  as  more  data  become  available. 
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Input  Requirements 

The  design  variables  required  for  the  new  design  method 
are  similar  to  those  required  by  the  current  FAA  procedure. 
These  variables  are  listed  below: 

a.  Traffic  in  aircraft  departure  per  year 

b.  Aircraft  type 

c.  Gross  aircraft  load,  P 

d.  Modulus  of  rupture  of  the  concrete,  Mr  psi 

e.  Modulus  of  subgrade  reaction,  K  in  pci 

The  strength  property  of  the  CTB  (i.e.  modulus  of 
rupture)  is  an  additional  input  required  by  the  new  method. 

The  factorial  design  performed  using  the  BISAR  computer 
program  required  elastic  material  properties  to  compute  the 
tensile  stresses  in  both  rigid  layers.  Therefore, 
relationships  were  needed  to  transform  K  of  the  subgrade  and 
Mr  of  the  CTB  into  elastic  moduli.  These  relationships  are 
listed  below: 

K  *  10(log[Esub]-1.415)/1.284  (6) 


Ectb  =  [Ln  (UC)  -  7.0]  /  3.597E-4  , (Ksi)  (7) 

Mrctb  -  218.2  EXP^3*597E”07*Ectb)  (8) 

Esub  =  subgrade  resilient  modulus,  psi 

UC  =  unconfined  compressive  strength  of  CTB,  psi 

Ectb>  ^ctb  =  modulus  of  elasticity  and  modulus  of  rupture 
of  the  CTB,  psi 

K  =  Modulus  of  subgrade  reaction,  pci 

An  example  design  is  performed  for  the  following 
factors: 


a.  Gear  type  =  dual  tandem 

b.  MR  of  concrete  =  700  psi 


c.  K  of  subgrade  =  150  pci 

d.  CTB  compressive  strength  =  2,000  psi 
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e.  Number  of  annual  departures  «  25,000 

f .  Gross  aircraft  load  »  380  kips 

The  modulus  of  rupture  of  the  CTB  is  related  to  the 
compressive  strength  by  the  following  relationship: 

MRctb  *  UC/5  (9) 

The  next  step  is  to  determine  the  concrete  thickness 
using  FAA  rigid  pavement  design  procedure.  The  procedure 
for  rigid  pavement  design  is  presented  in  an  Advisory 
Circular  (AC  No.  150/5320-6C)  dated  December  7,  1978.  Using 
the  input  requirement  stated  above,  with  the  exception  of 
the  CTB  strength  property,  a  concrete  thickness  can  be 
selected.  A  design  chart  is  used  for  each  gear  type. 

Knowing  the  design  modulus  of  rupture  of  the  concrete,  the 
gross  load  of  the  aircraft,  the  modulus  of  subgrade 
reaction,  and  the  number  of  annual  departures,  the  design 
thickness  can  be  determined  from  the  selected  design  chart 
(see  example  in  Figure  41) . 

Evaluation  of  the  tensile  stress  condition  in  the 
concrete  slab  of  a  one-rigid  layer  pavement  system  is  the 
next  step.  After  selecting  the  slab  thickness,  the  maximum 
tensile  stress  can  be  determined  using  the  regression 
equations  developed  in  this  study.  Three  equations  are 
provided  for  three  different  gear  types:  single  wheel,  dual 
wheel,  and  dual  tandem  gears.  Several  computational  steps 
are  required  prior  to  entering  the  regression  equations 
since  some  of  the  variables  in  these  equations  are  functions 
of  the  input  variables  described  in  the  previous  paragraph. 
Such  functions  include  the  radius  of  relative  stiffness,  L 

a.  Compute  L 
L  =  4  Econc*Kconc3 

12 (1-u2)K 

where 


ECOnc  =  Elastic  modulus  of  the  concrete,  psi 
HCOnc  =  Design  concrete  thickness,  in 

u  =  concrete  poisson's  ration  (=  0.15) 

K  =  modulus  of  subgrade  reaction,  pci 
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Figure  41.  Typical  design  chart  used  to  select  the  required  concrete  slab  thickness  for 

given  loading  and  soil  support  conditions  (FAA  Des.  Meth.,  AC  No.  150/5320-6C) 


b.  Use  the  following  relationships  to  determined  subgrade 
properties  K  -  io(lo9tEsubl  ”  1*415) /l. 284  ,  or 

Esub  -  lO*1*284  lo9CK)  +  1 .415) 

SU.  Determine  maximum  tensile  stress  in  one  rigid  layer 
system  using  regression  equations 

L  =  63.8  in 

Esub  =  16,185  psi 

Sc  =  325  psi  (from  regression  equations) 

d.  Assume  concrete  and  CTB  thicknesses  and  evaluate  stress 
condition  in  the  concrete  slab  of  the  two-rigid  layer 
pavement  system 

Initial  rigid  layer  thicknesses  are  selected  and  used 
to  compute  the  maximum  tensile  stress  in  the  concrete  slab. 
Equation  2  is  used  in  the  process.  Several  intermediate 
computational  steps  are  required  prior  to  entering  the 
regression  equations.  In  addition,  the  gear  type  is 
required  to  select  the  appropriate  regression  coefficients. 
Some  of  the  computational  steps  involved  are  listed  below: 

(1) .  Compute  modular  ratios: 


Eri  -  Econc/Ectb 

Er2  -  Ec^/Esub 

(2) .  Compute  radius  of  relative  stiffness  for  each  of  the 
rigid  layers  assuming  each  layer  is  supported  directly 
by  the  subgrade: 

Ll=  4  Econc  *  Hconc3 
12 ( 1-u2) K 

l2  =  Ectb  *  Hctb 

12 (1-u2) K 

where , 

Econc'  Ectb  =  modulus  of  elasticity  of  the  concrete 
and  CTB  in  psi,  respectively 
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u  -  Poisson's  ratio  of  the  concrete  and  CTB  materials. 

-  1, 684 , 000  psi 
Assume  H^-  11.5  -  10 

Er^  «  2.394 
Er2  -  103.1 
Lj  ®  43.1  in. 

1*2  *  31.3  in. 

sconc  88  249  Psi  sctb  "  192  Psi 

e.  Check  for  equality  of  maximum  tensile  stress  between  the 
rigid  pavement  systems 

The  maximum  tensile  stress  in  the  two  pavement 
structures  are  compared.  If  these  stresses  are  within 
+2.5  psi,  then  the  selected  concrete  and  CTB  thicknesses 
will  provide  equivalent  performance  to  the  one-rigid  layer 
pavement  system  designed  using  the  FAA  procedure.  If  the 
stress  difference  is  outside  the  allowable  range,  new 
concrete  and  CTB  thicknesses  are  assumed  and  step  d  is 
repeated.  This  process  is  continued  until  convergence  is 
obtained.  In  general,  convergence  is  achieved  within  three 
iterations.  The  result  of  this  step  is  the  selection  of 
layer  thicknesses  meeting  the  stress  equality  criteria.  An 
example  of  such  criteria  is  as  follows: 


FAA  design:  Hcom;  *  19.4  in. 

Stress  in  concrete  layer  =  235  psi 
Select  :  Hconc  -11.5  in. 

Hctjj  from  new  design  =  10.0  in. 

Stress  in  concrete  layer  =  249  psi 

249  >235  psi  =>  repeat  step  d 
Select:  Hconc  =  11.4  in,  Hctb  =  11.0  in. 

Stress  in  concrete  =  234  <  235  psi  =>  acceptable, 
f.  Check  for  maximum  tensile  stress  in  CTB 
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The  new  design  method  is  based  on  a  fixed  condition  of 
the  CTB.  It  is  assumed  that  the  condition  of  the  CTB  will 
remain  constant  over  the  life  of  the  pavement.  Therefore,  a 
limit  is  set  against  fatigue  failure  of  the  CTB.  For  the 
lack  of  any  other  published  data,  the  PCA  fatigue  failure 
for  concrete  pavement  is  criterion  applied  to  the  CTB.  This 
criterion  states  that  the  slab  can  resist  an  infinite  number 
of  load  applications  provided  the  stress  ratio,  i.e.  maximum 
tensile  stress  in  the  CTB  divided  by  the  modulus  of  rupture 
of  the  CTB,  remains  below  0.5. 

The  following  failure  criterion  obtained  from 
Thompson10  is  used  to  determine  the  number  of  allowable  load 
repetitions  if  the  stress  ratio  is  greater  than  0.5: 

Log  (N)  *  [(0.972  -  SR) /0 . 0825 ]  (10) 

where 


N  =  allowable  load  repetitions 

SR  =  stress  ratio,  stress  divided  by  MR 

It  is  recommended  that  the  stress  ratio  be  always  kept 
under  0.5  for  the  final  selection  of  the  layer  thicknesses. 
At  the  present  time,  the  limit  on  the  value  of  the  stress 
ratio  is  not  a  limiting  criterion.  However,  as  more  data 
becomes  available  regarding  the  fatigue  behavior  of  the  CTB, 
the  limiting  criterion  should  become  a  controlling  factor  in 
the  final  selection  of  the  layer  thicknesses.  This  is  shown 
in  the  following  example. 


Hconc  —  11*5  in. ,  —  11.0  m. 

sctb  =  181.2  psi 

SR  =  .453  =>  acceptable. 


Disadvantages  of  the  new  method 

Two  disadvantages  were  noted  in  the  application  of  the 
new  design  method.  The  first  disadvantage  is  the  tedious 
calculations  involved  in  the  manual  solution  of  the 
regression  equations.  This  disadvantage  stems  from  the 
large  number  of  terms  included  in  the  regression  equations. 
In  addition,  the  trial-and-error  procedure  involved  in  the 
final  selection  of  the  layer  thicknesses  in  order  to  meet 
the  allowable  stress  difference  between  the  one-and  two- 
rigid  layer  pavement  systems  is  very  time-consuming. 
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To  eliminate  this  disadvantage,  microcomputer  programs 
have  been  developed.  A  description  of  the  programs  is 
provided  in  the  next  section. 

The  second  disadvantage  is  that  of  layer  thickness 
selection.  The  selection  of  the  layer  thicknesses  is  not 
unique  since  the  stress  equality  between  the  two  pavement 
systems  can  be  obtained  for  different  layer  thicknesses  and 
layer  material  properties.  Therefore,  the  optimum  solution 
will  depend  on  the  judgment  of  the  engineer  performing  the 
design. 

Computer  Programs 

The  time  required  in  the  manual  selection  of  the 
concrete  slab  and  CTB  thicknesses  using  the  new  design 
procedure  is  examined.  This  fact  has  been  recognized  in  the 
previous  paragraph  as  one  of  the  main  disadvantages  of  the 
new  design  methods.  To  alleviate  this  problem,  two  computer 
codes  are  provided  to  automate  the  new  design 

A  description  of  both  programs  is  provided  in  the 
following  sections. 

General  Description 

The  two  computer  codes,  CTBDES  and  CTBEVL ,  are 
interactive  programs  which  are  used  to  design  or  evaluate 
two-rigid  layer  pavement  systems.  The  programs  solve  the 
regression  equations  obtained  from  the  large  factorial 
design.  The  source  codes  are  written  in  the  BASIC  computer 
language  and  are  suitable  for  use  with  an  IBM-compatible 
microcomputer.  Both  source  codes  are  included  in  Appendix 
C.  Minor  changes  may  be  required  to  the  source  codes, 
particularly  in  input/output  operations  such  as  opening  and 
closing  files,  in  order  to  be  used  on  other  computer 
systems . 

DESCRIPTION  OF  CTBDES 

CTBDES  is  used  to  select  the  concrete  and  CTB  slab 
thicknesses.  The  main  difference  between  CTBDES  and  CTBEVL 
is  that  the  first  program  provides  all  concrete-CTB 
thickness  combinations  that  satisfy  the  maximum  tensile 
stress  equality  criteria  described  earlier  in  this  report. 
All  concrete  and  CTB  thicknesses  within  the  valid  range  of 
the  regression  equations  are  evaluated.  Figure  42  presents 
the  flow  chart  of  CTBDES.  The  same  steps  involved  in  the 
hand  computations,  which  were  described  earlier  in  this 
SECTION,  are  followed  by  the  program. 

The  input  requirements  include  the  concrete  slab 
thickness  obtained  from  the  current  FAA  design  procedure  for 
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a  given  aircraft  type  and  design  gross  load,  and  the 
subgrade,  and  CTB  strength  properties.  Figure  43 
illustrates  an  example  of  the  input  requirements  for  CTBDES . 
The  computer  code  has  been  structured  such  that  more  than 
one  strength  property  can  be  used  to  characterize  a  given 
layer  of  the  pavement  system.  For  example,  the  modulus  of 
subgrade  reaction,  K,  or  the  resilient  modulus  of  the 
subgrade,  2SUb,  could  be  entered.  Similarly,  the 
compressive  strength  or  the  elastic  modulus  of  the  CTB  could 
be  entered. 

Figure  44  presents  an  example  of  the  output  obtained 
from  CTBDES.  The  user  input  is  echoed  in  the  top  portion  of 
the  output.  The  subgrade  and  base  properties  which  were  not 
entered  by  the  user  and  were  needed  for  the  solution  of  the 
regression  equations  are  computed  internally  by  the  program 
using  the  relationships  presented  earlier  in  this  section. 

The  concrete  and  CTB  thickness  combinations  are 
presented  in  the  lower  portion  of  the  output.  The  concrete 
slab  thickness  is  increased  by  1/2  in.  for  each  step,  and 
the  required  CTB  thickness  is  determined.  In  addition,  the 
maximum  tensile  stress  in  the  concrete  slab  of  the  rigid 
layer  pavement  system  (predicted)  is  compared  with  that  of 
the  one  rigid  layer  pavement  system  (design) . 

The  stress  ratio  in  the  CTB  is  also  printed  along  with 
the  resulting  allowable  load  repetitions  (i.e.  passes). 

DESCRIPTION  OF  CTBEVL 

Figure  45  presents  the  flow  chart  followed  by  CTBEVL. 
The  main  purpose  of  CTBEVL  is  to  select  the  required 
concrete  thickness  for  an  existing  pre-selected  CTB 
thickness.  CTBDES  was  restructured  as  indicated  on  the  flow 
charts  to  achieve  this  purpose.  In  general,  similar  outputs 
are  obtained  by  both  programs,  except  that  only  one 
concrete-CTB  thickness  combination  exists  for  CTBEVL. 

Design  Comparison 


The  new  design  method  developed  to  rationally  select 
the  required  concrete  and  CTB  thicknesses  was  presented  in 
Section  5.  This  section  describes  a  comparison  between  the 
results  obtained  by  the  new  design  method  and  those 
determined  using  the  current  FAA  design  method.  The  purpose 
of  this  comparison  is  to  validate  the  results  obtained  by 
the  new  method  and  to  study  its  sensitivity  to  changes  in 
the  design  parameters. 


WHICH  SUBGRADE  PROPERTY  DO  YOU  WANT  TO  USE  (K  OR  E>  ~  K 
ENTER  MODULUS  OF  SUBGRADE  REACTION,  K  <PC:>~  ISO 
WHICH  CTB  MATERIAL  PROPERTY  DO  YOU  WANT  TO  USE  <FC  OR  E> FC 
ENTER  COMPRESSIVE  STRENGTH  OF  CTB  <P SI)"  COCO 

ENTER  DESIGN  CONCRETE  THICKNESS  OBTAINED  FROM  FAA  PROCEDURE  (IN. )  ?  19.4 

ENTER  GEAR  TYPE  (SW.DT.DW)'?  DT 

ENTER  AIRCRAFT  GROSS  LOAD,  P  (KIPS)?  "80 


Figure  43.  Example  of  input  requirements  for  the  CTBDES  program 


5UBGF 

;AGE  PROPERTIES 

-  ESUB 

(  K.S  I )  = 

16. 18457 

K 

.(PCI)  = 

150 

BASE 

PROPERTIES  - 

ECTB 

(KSI)  = 

1694.584 

FC 

(PSI)  = 

2000 

CONCRETE  THICKNESS  AS  CALCULATED  BY  FAA  DESIGN,  HC  (IN)  =  19.40 

GEAR  TYPE  :  DT 

GROSS  LOAD  OF  DESIGN  AIRCRAFT,  P  (KIPS)  =  380 


DETERMINATION  OF  CONCRETE  AND  CTB  THICKNESS  COMBINATIONS 


THICKNESS  (IN) 

STRESS 

(PSI) 

STRESS 

ALLOWABLE 

CONC. 

CTB 

PREDICTED 

DESIGN 

RAT  ID 

PASSES 

8.0 

12.9 

237 

235 

0.727 

9.27E+02 

8.5 

12.8 

237 

235 

0.  673 

4. 17E+03 

9.0 

12.6 

237 

235 

0.626 

1.57E+04 

9.5 

12.5 

235 

235 

0.584 

5. 1 1E+04 

10.0 

12.2 

234 

235 

0.546 

1 . 47E-*-05 

10.5 

11.9 

234 

235 

0.512 

3.80E+05 

11.0 

11.2 

237 

235 

0.480 

9. 1 1E+05 

11.5 

11.0 

234 

235 

0.453 

1 . 97E+06 

12.0 

10.6 

233 

235 

0.427 

4.00E+06 

12.5 

10.0 

233 

235 

0.404 

7.  69E+06 

13.0 

9.3 

235 

235 

0.382 

1 . 42E+07 

13.5 

8.8 

233 

235 

0.362 

2.45E+07 

14.0 

8.0 

235 

235 

0.344 

4. 14E+07 

14.5 

7.5 

233 

235 

0.327 

6.62E+07 

15.0 

6.9 

233 

235 

0.311 

1 . 03E+08 

15.5 

6.2 

233 

235 

0.296 

1 . 57E+08 

16.0 

5.5 

233 

235 

0.282 

2.34E+08 

16.5 

4.7 

233 

235 

0.268 

3. 41E+08 

FOR  THE  LAST  ITERATION 

HCTB 

IS  <  4  IN.  -  USE 

HCTB=  4  IN. 

AND  HC=  17 

IN. 

Figure 

44.  Example  of 

the  output 

obtained 

using  CT 

BDES  program 

Figure  45.  Flow  chart  for  the  CTBEVL  program 


65 


Approach 

A  factorial  design  was  selected  to  perform  the 
comparison  of  the  required  concrete  slab  thickness  as 
obtained  by  the  two  design  methods.  The  factoxs  and  the 
levels  of  factors  included  in  the  factorial  design  are 
listed  below: 

a.  Gear  type:  dual  gear  and  dual  tandem 

b.  CTB  thickness:  6,  8  and  10  in. 

c.  CTB  elastic  modulus:  250,  500,  and  1,000  ksi 

d.  Subgrade  modulus:  5,  10,  20,  and  30  ksi 

e.  Traffic  level:  3,000,  15,000,  and  25,000  annual 
departures 

The  current  FAA  design  method  cannot  account  for 
changes  in  the  strength  characteristics  of  the  CTB  in  its 
determination  of  the  subgrade  support  value.  For  this 
reason,  a  third  design  method  was  added  to  the  comparison  in 
order  to  evaluate  better  the  capabilities  of  the  new  design 
method.  The  third  design  method  uses  a  relationship 
developed  by  the  US  Army  Corps  of  Engineers  (CE)  to  compute 
a  composite  K  of  the  subgrade  and  CTB  which  will  then  be 
used  in  the  current  FAA  design  method.  The  CE  K-composite 
charts  adjust  for  different  CTB  strength  properties  and 
thicknesses.  However,  K-composite  cannot  exceed  500  pci 
regardless  of  the  strength  or  thickness  of  the  CTB. 

A  total  of  504  designs  was  performed  as  a  part  of  this 
comparison. 

New  Design  Method 

The  computer  code  CTBEVL  was  used  to  perform  the 
computations.  The  results  are  listed  in  Table  6.  These 
results  indicate  that  the  new  method  is  sensitive  to  changes 
in  the  levels  of  the  factors  considered  in  the  design.  The 
required  concrete  slab  thickness  increases  as  the  resilient 
modulus  of  the  subgrade  decreases  and  the  traffic  level 
increases.  This  relationship  holds  true  for  all 
combinations  of  the  CTB  thickness  and  CTB  elastic  modulus  as 
well  as  for  the  two  gear  types.  In  addition,  for  a  given 
gear  type,  traffic  level,  and  subgrade  strength,  the 
required  concrete  slab  thickness  decreases  as  the  CTB 
thickness  and  CTB  elastic  modulus  increase. 
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The  new  design  method  does  not  set  any  limits  on  the 
structural  value  of  the  CTB  provided  that  the  CTB  and 
subgrade  material  properties,  the  concrete,  and  CTB 
thickness  combinations  are  within  the  valid  range  of  the 
regression  equations.  The  selected  concrete  slab 
thicknesses  were  based  on  providing  a  CTB  stress  ratio  of 
0.5  or  less.  Therefore,  fatigue  failure  is  not  anticipated 
to  occur  in  the  CTB  during  the  design  life  of  the  pavement 
structure . 

FAA  Method  (K-composite  determined  from  FAA  chart) 

This  design  method  was  performed  according  to  the 
procedure  included  in  the  Advisory  Circular  AC  150/5320-6C, 
dated  December  7,  1978.  Using  the  K-composite  chart 
provided  in  this  manual,  it  is  obvious  that  the  design 
method  is  insensitive  to  changes  in  the  elastic  modulus  of 
the  CTB.  However,  all  other  relationships  described  in  the 
previous  paragraph  also  apply  to  this  method.  The  results 
of  this  analysis  are  presented  in  Table  7. 

FAA  Method  (K-composite  determined  from  CE  chart) 

The  K-composite  chart  developed  by  the  CE  was  used  in 
this  method  along  with  the  design  chart  provided  in 
AC  150/5320-6C  for  each  aircraft  type  in  order  to  determine 
the  required  concrete  slab  thicknesses.  The  results 
obtained  from  this  method  are  listed  in  Table  8. 

This  method  is  insensitive  to  changes  in  subgrade  and 
CTB  material  properties  for  thicknesses  and  moduli 
combinations  producing  K-composite  values  greater  than 
500  pci.  The  occurrence  of  this  constraint  in  actual  design 
cases  is  illustrated  by  the  results  of  the  analysis  for  CTB 
thicknesses  of  8  and  10  in.,  subgrade  moduli  of  20  and 
30  ksi,  and  CTB  moduli  of  500  and  1,000  ksi.  For  a  given 
traffic  level,  tne  required  concrete  slab  thickness  is  the 
same  for  any  combination  of  the  three  variables  listed 
above.  Since  the  CE  method  is  not  dependent  upon  gear  type, 
this  condition  is  exhibited  for  both  gears. 
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Table  6 


Required  Concrete  Thickness  as  Determined  by  Hew  Deaton  Method 


Duel  Beir  Dull  Tindee 

Subgnde  Elutic  Modulus 
(kn) 


CTB  CTB  Traffic 

Thick.  Modulus  Uvei 

(in.)  (ksi)  (Dct./Yr.) 


3000 

250  15000 

25000 

3000 

6  500  15000 

25000 

3000 

1000  l”  00 

25000 

3000 

250  15000 

25000 

3000 

8  500  15000 

25000 

3000 

1000  15000 

25000 

3000 

250  15000 

25000 

3000 

10  500  15000 

25000 

3000 

1000  15000 

25000 


_5_ 

10 

JL 

30 

15.8 

14.6 

13.7 

13.4 

17.4 

16.2 

15.2 

14.7 

18.1 

16.7 

15.7 

15.3 

14.8 

13.8 

’2.9 

12.7 

16.5 

15.3 

14.4 

14.1 

17.1 

15.8 

14.9 

14.6 

13.9 

12.9 

12.2 

12.1 

15.6 

14.4 

13.7 

13.5 

16.2 

14.9 

14.2 

14 

15.3 

14.1 

13.7 

13.1 

16.9 

15.7 

14.8 

14.5 

17.5 

16.2 

15.3 

15 

14 

13.1 

12.3 

12.2 

15.7 

14.6 

13.8 

13.7 

16.3 

15.1 

14.3 

14.1 

12.8 

11.9 

11.3 

11.3 

14.5 

13.5 

12.8 

12.7 

15.2 

14 

13.4 

18.2 

14.6 

13.6 

12.8 

12.7 

16.3 

15.2 

14.4 

14.1 

16.9 

15.7 

14.8 

14.6 

13.3 

12.3 

11.7 

11.7 

14.9 

13.9 

13.2 

13.1 

15.6 

14.4 

13.8 

13.6 

11.7 

10.8 

10.3 

10.4 

13.4 

12.4 

11.9 

11.9 

14 

13 

12.4 

12.3 

5  10  20  30 


16.6 

15.2 

14.1 

13.3 

18.4 

16.9 

15.7 

14.7 

19.2 

17.4 

16.1 

15.2 

15.4 

14.1 

13.2 

12.4 

17.1 

15.8 

14.6 

13.9 

17.9 

16.2 

15.2 

14.2 

14.3 

13 

12.2 

11.6 

16 

14.7 

13.8 

13 

16.7 

15.3 

14.2 

13.4 

15.9 

14.5 

13.7 

12.9 

17.6 

16.2 

15.3 

14.4 

18.4 

16.7 

15.7 

14.7 

14.4 

13.2 

12.4 

11.8 

16.1 

14.9 

14 

13.3 

16.8 

15.4 

14.4 

13.7 

12.8 

11. B 

11.2 

10.5 

14.6 

13.5 

12.6 

12 

15.4 

14 

13.2 

to  3 

15.2 

13.9 

13.1 

12.* 

16.9 

15.6 

14.6 

13.9 

17.6 

16.1 

15.2 

14.3 

13.4 

12.2 

11.6 

11.1 

15.2 

14 

13.2 

12.4 

15.9 

14.5 

13.7 

12.8 

11.4 

10.4 

9.9 

9.4 

13.2 

12.1 

11.5 

10.8 

13.9 

12.6 

12 

11.3 
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Table  7 


Required  Concrett  Thickness  Based  on  E-Composite  as  Determined 


bv  New  Design  Hcthod 


CTB 

Thick. 

(in.) 


6 


10 


CTB  Traffic  Dual  Bear 


Dual  Tandee 


Modulus  Level 
Iksi)  (Dep./Yr.) 


Subgrade  Elastic  Modulus 
(ksi) 


_5_ 

JO 

20 

30 

5 

10 

20 

30 

3000 

15.1 

14.3 

13.6 

13.4 

15.7 

14.1 

13.1 

12.7 

250 

15000 

16.6 

15.7 

15.1 

14.8 

17.2 

15.4 

14.4 

14 

25000 

17.1 

16.1 

15.5 

15.2 

17.8 

15.9 

14.8 

14.4 

3000 

15.1 

14.3 

13.6 

13.4 

15.7 

14.1 

13.1 

12.7 

500 

15000 

16.6 

15.7 

15.1 

14.8 

17.2 

15.4 

14.4 

14 

25000 

17.1 

16.1 

15.5 

15.2 

17.8 

15.9 

14.8 

14.4 

3000 

15.1 

14.3 

13.6 

13.4 

15.7 

14.1 

13.1 

12.7 

1000 

15000 

16.6 

15.7 

15.1 

14. B 

17.2 

15.4 

14.4 

14 

25000 

17.1 

16.1 

15.5 

15.2 

17. B 

15.9 

14.8 

14.4 

3000 

14.7 

14 

13.4 

13.2 

14. B 

13.5 

12.7 

12.2 

250 

15000 

16.1 

15.4 

14.B 

14.6 

16.3 

14.8 

14 

13.4 

25000 

16.6 

15.9 

15.2 

15 

16.8 

15.2 

14.4 

13.8 

3000 

14.7 

14 

13.4 

13.2 

14.8 

13.5 

12.7 

12.2 

500 

15000 

16.1 

15.4 

14.  B 

14.6 

16.3 

14. B 

14 

13.4 

25000 

16.6 

15.9 

15.2 

15 

16.8 

15.2 

14.4 

13.8 

3000 

14.7 

14 

13.4 

13.2 

14.8 

13.5 

12.7 

12.2 

1000 

15000 

16.1 

15.4 

14.8 

14.6 

16.3 

14.8 

14 

13.4 

25000 

16.6 

15.9 

15.2 

15 

16.  B 

15.2 

14.4 

13. B 

3000 

14.2 

13.5 

13.2 

12.8 

14 

13 

12.5 

12 

250 

15000 

15.5 

15 

14.5 

14,2 

15.3 

14.3 

13.8 

13.2 

25000 

16 

15.4 

14.9 

14.6 

15.8 

14.7 

14.2 

13.5 

3000 

14.2 

13.5 

13.2 

12.8 

14 

13 

12.5 

12 

500 

15000 

15.5 

15 

14.5 

14.2 

15.3 

14.3 

13.8 

13.2 

25000 

16 

15.4 

14.9 

14.6 

15.8 

14.7 

14.2 

13.5 

3000 

14.2 

13.5 

13.2 

12.8 

14 

13 

12.5 

12 

1000 

15000 

15.5 

15 

14.5 

14.2 

15.3 

14.3 

13.8 

13.2 

25000 

16 

15.4 

14.9 

14.6 

15.8 

14.7 

14.2 

13.5 
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Table  8 


Beeul  red  Concretf  Thickness. ftufd.egJ^gaePli^  >’  >4 nrined 
bv  Coroa  of  Engineers  Method 

CT8  CTB  Traffic  Dual  6ear  Du*l  Tandee 

Thick.  Modulus  level  Subg'ade  Elastic  Modulus 


lin.)  (ksi)  (Dep./Tr.) 

JL.  JO 


3000 

15.2 

14.4 

250 

15000 

16.8 

15.7 

25000 

17.2 

16.2 

3000 

14.9 

13.9 

500 

15000 

10.4 

15.3 

25000 

16.8 

15.7 

3000 

14.1 

13.2 

1000 

15000 

15.5 

14.6 

25000 

16 

15 

3000 

14.9 

14.1 

250 

15000 

16.4 

15.5 

25000 

16.8 

16 

3000 

14.7 

13.4 

500 

15000 

16.1 

14.8 

25000 

16.6 

15.2 

3000 

13.5 

13.2 

1000 

15000 

15.3 

14.6 

25000 

15.8 

15 

3000 

14.9 

14 

250 

15000 

16.4 

15.4 

25000 

16.8 

15.9 

3000 

14.7 

13.4 

500 

15000 

16.1 

14.8 

25000 

16.6 

15.2 

3000 

13.5 

12.5 

1000 

15000 

15 

13.8 

25000 

15.4 

14.2 

(ksii 


20 

30 

5 

i! 

J1 

30 

13.2 

12.5 

15.8 

14.4 

12.5 

11.3 

14.5 

13.8 

17.3 

15.8 

13.8 

12.4 

15 

14.2 

17.8 

16.3 

14.2 

12.8 

13.1 

12.5 

15.3 

13.7 

11.8 

11.3 

14.3 

13.8 

16.8 

15 

13 

12. J 

14.8 

14.2 

17.4 

15.5 

13.3 

12.7 

12.5 

12.5 

14 

12.1 

11.3 

11. J 

13.8 

13.8 

15.3 

13.4 

12.3 

12.3 

14.2 

14.2 

15.8 

13.8 

12.7 

12.7 

12.7 

12.5 

15.3 

14 

11.5 

11.3 

14 

13.8 

16.8 

15.3 

12.6 

12.3 

14.4 

14.2 

17.4 

15  B 

13 

12.7 

12.5 

12.5 

14.6 

12.7 

11.3 

11.3 

13.8 

13.6 

16.3 

14 

12.3 

12.3 

14.2 

14.2 

16.8 

14.4 

12.7 

12.7 

12.7 

12.5 

13.4 

Pe2 

11.3 

11.3 

14 

13.8 

14.8 

13.4 

12.3 

12.3 

14.4 

14.2 

15.2 

13.8 

12.7 

12.7 

12.5 

12.5 

15.3 

13.5 

11.5 

11.3 

13.8 

13.8 

16.8 

14.0 

32.6 

12.3 

14.2 

14.2 

17.4 

13.2 

13 

12.7 

12.5 

12.5 

14.8 

12.7 

11.3 

11.3 

13.8 

13.8 

16.3 

14 

12.3 

12.3 

14.2 

14.2 

16.8 

14.4 

12.7 

12.7 

12.5 

12.5 

13 

11.3 

11.3 

11.3 

13.8 

13.8 

14.3 

12.4 

12.4 

12.4 

14.2 

14.2 

14.7 

12.8 

12.6 

12.8 
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Difference  in  Results 


The  results  obtained  by  the  new  method  were  compared 
with  the  current  FAA  and  CE  methods  by  computing  the 
difference  between  the  required  concrete  slab  thicknesses 
(Hconc) new  ®in^s  (Hconc) current*  These  differences  are 
shown  in  Tables  10  and  11.  The  results  indicated  that  the 
current  FAA  design  method  tends  to  over  design  the  concrete 
slabs  for  the  cases  where  Ectb  is  equal  to  1,000  ksi.  On 
the  other  hand,  the  same  method  tends  to  underestimate  the 
required  slab  thickness  for  an  Ectb  of  250  ksi. 

In  general,  the  current  method  requires  a  greater 
concrete  thickness  than  the  newly  developed  method  in  64 
percent  of  the  cases  included  in  the  factorial  design.  The 
same  condition  occurred  when  the  concrete  thicknesses  were 
determined  using  the  K-composite  relationship  developed  by 
the  CE.  However,  only  62  percent  of  the  cases  included  in 
our  analysis  required  greater  concrete  thickness  than  that 
recommended  by  the  new  method. 

No  other  trends  were  observed  in  the  difference  between 
the  estimated  concrete  slab  thickness  using  the  three 
different  methods. 
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Table  9 


Dlfferer.ca  in  Thickness  of  Concrete  Slab  »8  Determined  bv  the  New  Oesiqn  Method 
and  the  C-CowootUe  wet  hod  (FAA  Hethodl 


CIS 

Thick. 

(in.) 


10 


CT8  TraHit 
Modulus  Level 
(rsi >  (Dep./Y r.) 


Duel  6eir  Dual  Tindei 

Subarade  Elastic  Modulus 

(irsi) 


JL 

11 

IL 

30 

x 

UL 

ii 

30 

3000 

0.7 

0.3 

0.1 

0 

0.9 

u 

i 

0.6 

250 

15000 

0.8 

0.5 

0.1 

-0.1 

1.2 

1.5 

1.3 

0.7 

25000 

1.1 

0.4 

0.2 

0.1 

1.4 

1.5 

1.3 

0.8 

3000 

*0.3 

-0.5 

-0.7 

-0.7 

-0.3 

0 

0.1 

-0.3 

500 

15000 

-0.1 

-0.4 

-0.7 

-0.7 

-0.1 

0.4 

0.2 

-0.1 

25000 

0.1 

-0.3 

-0.4 

-0.6 

0.1 

0.3 

0.4 

-0.2 

3000 

-1.2 

-1.4 

-1.4 

-1.3 

-1.4 

-1.1 

-0.9 

-1.1 

1000 

15000 

-t 

-1.3 

-1.4 

-1.3 

-1.2 

-0.7 

-0.6 

-1 

25000 

-0.8 

-1.2 

-1.3 

-1.2 

-1.1 

-0.6 

-0.6 

-1 

3000 

0.4 

0.1 

0.3 

-0.1 

1.1 

1 

1 

0.7 

250 

.5000 

0.8 

0.3 

0 

-0.1 

1.3 

1.4 

1.3 

1 

25000 

0.0 

0.3 

0.1 

0 

1.4 

1.5 

1.3 

0.9 

3000 

-0.7 

-0.9 

-1.1 

-1 

-0.4 

-0,3 

-0.3 

-0.4 

500 

15000 

-0.1 

-0.8 

-l 

-0.9 

-0.2 

0.1 

0 

-0.1 

25000 

-0.3 

-0.8 

-0.9 

-0.9 

0 

0.2 

0 

-0.1 

3000 

-1.9 

-2.1 

-2.1 

-1.9 

-2 

-1.7 

-1.5 

-1.7 

1000 

15000 

-1.4 

-1.9 

-2 

-1.9 

-1.7 

-1.3 

-1.4 

-1.4 

25000 

-1.4 

-1.9 

-1.8 

3.2 

-1.4 

-1.2 

-1.2 

-1.5 

3000 

0.4 

0.1 

-0.4 

-0.1 

1.2 

0.9 

0.6 

0.4 

250 

15000 

0.8 

0.2 

-0.1 

-0.1 

1.4 

1.3 

0.8 

0.7 

25000 

0.9 

0.3 

-0.1 

0 

1.8 

1.4 

1 

0.6 

3000 

-0.9 

-1.2 

-1.5 

-1.1 

-0.4 

-0.8 

-0.9 

-0.9 

500 

15000 

-0.4 

-1.1 

-1.3 

-1.1 

-0.1 

-0.3 

-0.6 

-0.8 

25000 

-0.4 

-1 

-1.1 

-1 

0.) 

-0.2 

-0.5 

-0.7 

3000 

-2.5 

-2.7 

-2.9 

-2.4 

-2.6 

-2.6 

-2.6 

-2.6 

1000 

15000 

-2.1 

-2.4 

-2.4 

-2.3 

-2.1 

-2.2 

-2.3 

-2.4 

25000 

-2 

-2.4 

-2.5 

-2.3 

-1.9 

-2.1 

-2.2 

-2.2 
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Table  10 


CTB  CTB  Traffic  Dual  Gear  Duai  Tandee 

Thick.  Modulus  Level  Subgrade  Elastic  Modulus 

(in.)  (ksi)  (Oep./Yr.)  (ksi) 


J5_ 

20 

30 

_5_ 

i! 

20 

30_ 

3000 

0.6 

0.2 

0.5 

0.9 

0.8 

0.8 

1.6 

2 

250 

15000 

0.6 

0.5 

0.7 

0.9 

1.1 

1.1 

1.9 

2.3 

25000 

0.9 

0.5 

0.7 

1.1 

1.4 

1.1 

1.9 

2.4 

3000 

-0.1 

-0.1 

-0.2 

0.2 

0.1 

0.4 

1.4 

1.1 

500 

15000 

0.1 

0 

0.1 

0.3 

0.3 

0.8 

1.6 

1.6 

25000 

0.3 

0.1 

0.1 

0.4 

0.5 

0.7 

1.9 

1.5 

3000 

-0.2 

-0.3 

-0.3 

-0.4 

0.3 

0.9 

0.9 

0.3 

1000 

15000 

0.1 

-0.2 

-0.1 

-0.3 

0.7 

1.3 

1.5 

0.7 

25000 

0.2 

-0.1 

0 

-0.2 

0.9 

1.5 

1.5 

0.7 

3000 

0.* 

0 

1 

0.6 

0.6 

0.5 

2.2 

1.6 

250 

15000 

0.5 

0.2 

0.8 

0.7 

0.8 

0.9 

2.7 

2.1 

25000 

0.7 

0.2 

0.9 

0.8 

1 

0.9 

2.7 

2 

3000 

-0.7 

-0.3 

-0.2 

-0.3 

-0.4 

0.5 

1.1 

0.5 

500 

15000 

-0.4 

-0.2 

0 

-0.1 

-0.2 

0.9 

1.7 

1 

25000 

-0.3 

-0.1 

0.1 

-0.1 

0 

1 

1.7 

1 

3000 

-1.1 

-1.3 

-1.4 

-1.2 

-0.6 

-0.4 

-0.1 

-0.8 

1000 

15000 

-0.8 

-1.1 

-1.2 

-1.1 

-0.2 

0.1 

0.3 

-0.3 

25000 

-O.i 

-1 

-1 

4 

0.2 

0.2 

0.5 

-0.4 

3000 

-0.3 

-0.4 

0.3 

0.2 

-0.1 

0.4 

1.6 

1.1 

250 

15000 

-0.1 

-0.2 

0.6 

O.J 

0.1 

0.8 

2 

1.6 

25000 

0.1 

-0.2 

0.6 

0.4 

0.2 

0.9 

2.2 

1.6 

3000 

-1.4 

-1.1 

-0.8 

-0.8 

-1.4 

-0.5 

0.3 

-0.2 

500 

15000 

-1.2 

-0.9 

-0.6 

-0.7 

-l.l 

0 

0.9 

0.1 

25000 

-0.8 

-0.4 

-0.6 

-0.9 

0.1 

1 

0.1 

3000 

-1.8 

-1.7 

-2.2 

-2.1 

-1.6 

-0.9 

-1.4 

-1.9 

1000 

15000 

-1.6 

-  -1.4 

-1.9 

-1.9 

-1.1 

-0.3 

-0.9 

-1.6 

25000 

-1.4 

-1.2 

-1.8 

-1.9 

-0.8 

-0.2 

-0.8 

-1.5 
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CLOSURE 


Summary 

This  study  has  resulted  in  the  development  of  a  design 
procedure  for  the  determination  of  the  thickness  of  a  cement 
stabilized  base  under  a  rigid  pavement  system.  The 
procedure  has  been  developed  based  on  an  equivalent  interior 
stress  criteria.  Allowance  has  been  made  for  the 
determination  of  the  allowable  pass  level  on  the  CTB 
material. 

Two  computer  programs  have  been  developed  in 
conjunction  with  the  study.  The  first  program  CTBDES  is  the 
design  program  to  be  utilized  in  the  thickness  determination 
for  new  pavement  structures.  The  second  CTBEVAL  is  an 
evaluation  program  which  can  be  used  in  checking  designs 
performed  by  other  methods,  or  it  can  be  used  in  an 
evaluation  of  an  existing  airfield. 

The  procedures  have  been  developed  using  elastic  layer 
analysis  in  the  form  of  the  BISAR  program.  The  method  has 
been  structured  around  the  existing  FAA  design  procedures 
and  is  completely  compatible  with  the  current  design 
procedures.  The  comparison  of  the  new  design  procedures  to 
the  existing  methods  show  that  the  thickness  differences  are 
both  plus  and  minus.  However,  there  are  several  cases  in 
which  the  new  procedure  shows  a  substantial  reduction  in  the 
required  thickness  of  concrete.  These  cases  occur  primarily 
when  a  CTB  mixture  of  high  elastic  properties  are  utilized 
in  the  design  process.  The  opposite  is  also  true.  When  low 
materials  properties  are  used  the  new  design  requires  up  to 
2  in.  of  more  concrete.  This  is  a  significant  finding.  If 
the  present  design  methods  are  used  without  consideration 
for  the  strength  of  the  CTB  material,  substantial  under¬ 
designs  are  possible. 

From  a  design  point  of  view  the  benefit  to  this 
analysis  procedure  is  that  the  designer  is  allowed  to 
utilize  the  materials  properties  of  the  CTB.  In  areas  where 
CTB  is  a  common  construction  material  or  for  large  projects, 
the  reductions  in  concrete  thickness  and  replacement  with  a 
higher  quality  CTB  may  provide  an  opportunity  for  economic 
savings.  Also,  the  designer  is  provided  with  information 
regarding  the  allowable  pass  level  on  the  CTB  material. 

This  information  is  not  presently  provided  in  any 
conventional  design  method. 
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Recommendations  for  Further  study 

The  major  area  encountered  in  this  study  which 
requires  further  research  is  the  fatigue  properties  of  the 
CTB  materials.  The  literature  is  very  sparse  concerning 
this  aspect  of  the  material.  For  the  most  part  references 
are  made  to  compressive  or  tensile  strengths;  however ,  no 
major  studies  have  been  performed  on  the  fatigue  properties. 
There  is  ample  work  performed  on  lean  concrete,  but  the 
strength  levels  of  these  materials  are  generally  above  the 
CTB  materials  in  the  250  to  500  psi  level. 

It  is  recommended  that  a  study  of  the  CTB  fatigue 
properties  be  undertaken  to  further  augment  this  study.  As 
an  alternate,  the  pavements  studied  in  the  "high  traffic 
volume"  study  should  be  continually  monitored  for  their 
field  performance.  These  pavement  systems  are  Portland 
Cement  Concrete  (PCC)  over  CTB  systems  where  the  traffic  and 
materials  properties  have  been  defined.  The  field 
performance  of  these  pavements  could  be  incorporated  into 
this  process  through  the  use  of  the  CTBEVAL  program  and  a 
correlation  with  the  Pavement  Condition  Index  (PCI)  or 
linear  feet  of  cracking  parameter. 

This  study  has  shown  that  the  stress  values  in  a  three 
layer  system  are  predictable.  It  is  recommended  that  the 
methodology  be  extended  to  cover  the  design  of  other  base 
materials  such  as  bituminous  stabilized  materials. 
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APPENDIX  A 


RESULTS  OF  THE  STRESS  ANALYSIS 
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APPENDIX  B 


ILLUSTRATION  OF  RELATIONSHIPS  BETWEEN  THE  VARIABLES 
INCLUDED  IN  THE  FACTORIAL 
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Figure  Bl.  Relationship  between  He  and  stress  in  the  concrete 
slab  for  all  CIB  moduli 
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Figure  B2. 

Relationship  between  Hctb  and  stress 
slab  for  all  CIB  moduli 
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Figure  B5.  Relationship  between  Hctb  and  stress  in  the  concrete 
slab  for  all  subgrade  moduli 
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Figure  B6.  Relationship  between  Esub  and  stress  in  the 
concrete  slab  for  all  CTB  moduli 
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Figure  B7. 


Relationship  between  He  and  stress  in  the 
concrete  slab  for  all  subgrade  moduli 
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Figure  B8.  Relationship  between  Ectb  and  stress  in  the 
concrete  slab  for  all  subgrade  moduli 
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Figure  B9.  Relationship  between  Esub  and  stress  in  the  CTB 
for  all  concrete  slab  thicknesses 


Figure  BIO.  Relationship  between  Ectb  and  stress  in  the  CTB 
for  all  subgrade  moduli 
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Figure  B15.  Relationship  between  Hctb  and  stress  in  the  CTB 
for  all  subgrade  moduli 
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Figure  B16. 


Relationship  between  subgrade  moduli  ard  str 
in  the  CTB  for  all  CTB  moduli 
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Figure  B17. 

Relationship  between  concrete  slab  thickness 
and  stress  in  the  concrete  slab  for  all  Esub 
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Figure 

B18. 

Relationship  between  concrete  slab 
and  stress  in  the  concrete  slab  for 
evaluated  at  the  boundaries  of  the 

thickness 
all  Esub 
factorial 
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Figure  B19.  Relationship  between  CTB  thickness  and  stress  in 
the  CTB  for  all  subgrade  moduli 
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Figure  B20.  Relationship  between  CTB  thickness  and  stress  in 
the  CTB  for  all  subgrade  moduli  evaluated  at  the 
boundaries  of  the  factorial 
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Figure  B21.  Relationship  between  concrete  slab  thickness 
and  stress  in  the  concrete  slab  for  all  CTB 
moduli 
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Figure  D22.  Relationship  between  CTB  thickness,  and  strops  in 
the  concrete  slab  for  all  CTB  moduli 
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Figure  B25.  Relationship  between  CTB  thickness  and  stress  in 
the  concrete  slab  for  all  subgrade  moduli 
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Figure  026.  Relationship  between  subgradc  moduli  aid 
in  the  concrete  slab  for  all  CTB  moduli 
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Figure  B27.  Relationship  between  concrete  slab  thickness  and 

stress  in  the  concrete  slab  for  all  svibgrade  moduli 
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Figure  P28.  Relationship  betv/ecn  CTB  moduli  and  stress  in  the 
concrete  slab  for  all  subgrade  moduli 


Figure  B30.  Relationship  between  CTB  moduli  and  stress  in  the 
CTB  for  all  subgrade  moduli 


Dual  Wheel 
Esub  =  10  ksi 
Hctb  =  10  in. 


0  s  10  15 

i  i 

20  25 

He 

-a—  500  -*-iooo 

-1-2500 

Figure  B31.  Relationship  between  concrete  slab  thickness 
and  stress  in  the  CTB  for  all  CTB  moduli 
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Figure  B32.  Relationship  between  CTB  thickness  and  stress 
in  the  CTB  for  all  CTB  moduli 
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Figure  B35.  ^  Relationship  between  CTB  thickness  and  stress  in  the 
"  CIB  for  all  subgrade  moduli 


Figure  B36.  Relationship  between  suhgrade  moduli  and  stress 
in  the  CTB  for  all  CTB  moduli 
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Figure  B37.  Relationship  between  concrete  slab  thickness 
and  stress  in  the  concrete  slab  for  all  Esub 
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Figure  B38. 

Relationship  between  concrete  slab  thickness 
and  stress  in  the  concrete  slab  for  all  Esub 
evaluated  at  the  boundaries  of  the  factorial 
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Figure  B39.  Relationship  between  CTB  thickness 
in  the  CIB  for  all  subgrade  moduli 

and  stress 
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Relationship  between  CTB  thickness  and  stress 
in  the  CTB  for  all  subgrade  moduli  evaluated  at 
fhf»  boundaries  of  the  factorial 
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Figure  B41.  Relationship  between  concrete  slab  thickness 
and  stress  in  the  concrete  for  all  CTB  moduli 

Dual  Tandem 
Esub  =  10  ksi 
He  =  14  in. 


0  2  4  6  8  10  12  14  16  18  20  22  24 


Figure  B42.  Relationship  between  CTB  thickness  and  stress 
in  the  concrete  for  all  CTB  moduli 
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Figure  B45.  Relationship  between  CIB  thickness  and  stress 
in  the  concrete  for  all  subgrade  moduli 
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Figure  B46.  Relationship  between  subgrade  moduli  and  stress 
in  the  concrete  for  all  CTB  moduli 
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Figure  B47.  Relationship  between  concrete  thickness  and  stress 
in  the  concrete  for  all  subgrade  moduli 
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Figure  B48.  Relationship  between  CTB  moduli  and  stress  in 
the  concrete  for  all  subgrade  moduli 
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Figure  B49.  Relationship  between  concrete  thickness  and 
stress  in  the  CTB  for  all  subgrade  moduli 
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Figure  B50.  Relationship  between  CTB  moduli  and  stress 
in  the  CTB  for  all  subgrade  moduli 
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Figure  B52.  Relationship  between  CIB  thickness  and  stress 
in  the  CIB  for  all  CTB  moduli 
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Figure  B54.  Relationship  between  CTB  moduli  and  stress  in 
the  CTB  for  all  CTB  thicknesses 
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Figure  B55.  Relationship  between  CTB  thickness  and  stress  in 
the  CIB  for  all  subgrade  moduli 
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Figure  B56.  Relationship  between  subgrade  moduli  and  stress 
in  the  CTB  for  all  CTB  moduli 
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Figure  B57.  Relationship  between  concrete  thickness  and  stress 
in  the  concrete  slab  for  all  subgrade  moduli 
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Figure  B58.  Relationship  between  concrete  thickness  and  stress 
in  the  concrete  slab  for  all  suberade  moduli 
evaluated  at  the  boundaries  of  tne  factorial 
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Figure  B60.  Relationship  between  CTB  thickness  and  stress  in 
the  CTB  for  all  subgrade  moduli  evaluated  at  the 
boundaries  of  the  factorial 


APPENDIX  C 

SOURCE  CODES  FOR  CTBDES  AND  CTBEVL  PROGRAMS 


CTBDES  PROGRAM 


10  KEY  OFFtCLS 
20  LOCATE  10,5 

30  PRINT  "THIS  PROGRAM  PROVIDES  A  LIST  OF  THICKNESS  COMBINATIONS 
OF  CONCRETE" 

40  PRINT  "  SURFACE  AND  CEMENT  TREATED  BASE.  THE  SELECTED 
THICKNESS  COMBINATIONS" 

50  PRINT  "  ARE  BASED  ON  ACHIEVING  AN  EQUIVALENT'  MAXIMUM  TENSILE 
STRESS  AT  THE 

60  PRINT  "  BOTTOM  OF  THE  CONCRETE  SURFACE  LAYER.  STRESS  RATIOS 
IN  THE  BASE  ARE" 

70  PRINT  "  PROVIDED  FOR  EACH  COMBINATIONS.  A  STRESS  RATIO  OF 
0.5  OR  LESS  IS* 

80  PRINT  "  RECOMMENDED." 

90  LOCATE  24,25 

100  PRINT  "PRESS  ANY  KEY  TO  CONTINUE" 

110  A$«INKEY$ 

120  IF  A$»""  THEN  GOTO  110  ELSE  130 

130  CLS 

140  LOCATE  5 

150  PRINT  "INPUT  REQUIREMENTS  INCLUDE  SUBGRADE  PROPERTY,  THICKNESS 
OF  CONCRETE" 

160  PRINT  "SURFACE  LAYER  AS  DETERMINED  BY  FAA  DESIGN  PROCEDURE, 
LOAD  CONDITION," 

170  PRINT  "AND  THE  CEMENT  TREATED  BASE  MATERIAL  PROPERTY.  THE 
STRESSES  USED" 

180  PRINT  "TO  ACHIEVE  THE  EQUIVALENCY  BETWEEN  THE  TWO  PAVEMENT 
SYSTEM  ARE" 

190  PRINT  "COMPUTED  USING  LINEAR  REGRESSION  EQUATIONS.  THE  VALID 
RANGE  OF" 

200  PRINT  "VARIABLES  LISTED  BELOW.  NO  ERROR  MESSAGES  WILL  BE 
PRINTED  IF" 

210  PRINT  "THE  SELECTED  VARIABLE  VALUE  IS  OUTSIDE  THE  VALID  RANGE 
OF  THE  EQUATIONS." 

220  LOCATE  13 

230  PRINT  "  SUBGRAGE  RESILIENT  MODULUS  =  5  TO  30  KSI" 

240  PRINT  "  MODULUS  OF  SUBGRADE  REACTION  *  60  TO  240  PCI" 

250  PRINT  "  THICKNESS  OF  CONCRETE  LAYER  *  8  TO  26  INCHES" 

260  PRINT  "  THICKNESS  OF  CTB  -  4  TO  22  INCHES" 

270  PRINT  "  LOAD  CONDITION  -SINGLE  WHEEL  GEAR,  DUAL  GEAR," 

280  PRINT  "  AND  DUAL  TANDEM  GEAR" 

290  PRINT  "  CTB  RESILIENT  MODULUS  «  250  TO  2,500  KSI" 

300  PRINT  "  CTB  COMPRESSIVE  STRENGTH  -  1,200  TO  2,700  PSI" 

310  LOCATE  24,25 

320  PRINT  "PRESS  ANY  KEY  TO  CONTINUE" 

330  A$-INKEY$ 

340  IF  A$-""  THEN  GOTO  330  ELSE  350 
350  CLS 

360  OPTION  BASE  1:DIM  A(3,7),  B(3,13),  C(3,13) 

370  GOSUB  860 
380  LOCATE  10 


cs 


390  INPUT" DO  YOU  WANT  TO  SEND  OUTPUT  TO  SCREEN.  PRINTER,  OR  FILE 
(S,  P,OR  F) " ;C$ 

400  IF  C$-"S"  OR  C$="s"  THEN  430 

410  IF  C$-"f"  OR  C$«"F"  THEN  GOTO  440 

420  IF  C$«"p"  OR  C$*"P"  THEN  GOTO  460  ELSE  390 

430  OPEN  "SCRN: "  FOR  OUTPUT  AS  #l:GOTO  470 

440  INPUT"NAME  OF  YOUR  OUTPUT  FILE";STT$ 

450  OPEN  STT$  FOR  OUTPUT  AS  #l:GOTO  470 
460  OPEN  "LPT1"  FOR  OUTPUT  AS  #1 
470  CLS 
480  LOCATE  10 

490  INPUT"WHICH  SUBGRADE  PROPERTY  DO  YOU  WANT  TO  USE  (K  OR  E)";A$ 
500  IF  A$*"K"  OR  A$="k"  THEN  GOTO  540 

510  IF  A$«"E"  OR  A$«"e"  THEN  GOTO  520  ELSE  GOTO  490 

520  INPUT"ENTER  SUBGRADE  RESILIENT  MODULUS,  E 
(KSI) " ;ESUB:BSUB-ESUB*1000 

530  K**lOA(((LOG(:i!SUB)/LOG(10))-1.415)/1.284)  :  GOTO  560 
540  INPUT"ENTER  MODULUS  OF  SUBGRADE  REACTION,  K  (PCI)";K 
550  ESUB*10A (1.284* ( LOG (K) /LOG (10) ) +1.415) 

560  INPUT "WHICH  CTB  MATERIAL  PROPERTY  DO  YOU  WANT  TO  USE  (FC  OR 
E) " ;B$ 

570  IF  B$-"FC"  OR  B$«"fc"  THEN  GOTO  610 
580  IF  B$-"E"  OR  B$-"e"  THEN  GOTO  590  ELSE  GOTO  560 
590  INPUT" ENTER  MODULUS  OF  CTB  (KSI) " ;ECTB:ECTB»ECTB*1000 
600  UC*1091 . 2*EXP ( 3 . 596552E-07 *ECTB) : GOTO  630 
610  INPUT "ENTER  COMPRESSIVE  STRENGTH  OF  CTB  (PSI)";UC 
620  ECTB- (LOG (UC) -LOG (1091 . 2 ) ) /3 . 596552E-07 
630  INPUT" ENTER  DESIGN  CONCRETE  THICKNESS  OBTAINED  FROM  FAA 
PROCEDURE  (IN.)";HCD 
640  ER-40000001/ESUB 

650  L-(4000000l*(HCDA3)/(11.73*K) ) A.25 
660  INPUT"ENTER  GEAR  TYPE  (SW, DT, DW) " ;GT$ 

665  INPUT"ENTER  AIRCRAFT  GROSS  LOAD,  P  (KIPS)  "  ;P:  P=»P*1000 
670  CLS 

680  IF  C$«"p"  OR  C$«"P"  OR  C$«"F"  OR  C$="f"  THEN  GOTO  690  ELSE 
GOTO  710 
690  LOCATE  12,15 

700  PRINT  "PLEASE  WAIT  FOR  OUTPUT  FILE  OR  PRINTER" 

710  PRINT  #1, "SUBGRAGE  PROPERTIES  -  ESUB  (KSI)  -  ",ESUB/1000 

720  PRINT  #1,"  K  (PCI)  *  ",K 

730  PRINT  #1, "BASE  PROPERTIES  -  ECTB  (KSI)  -",ECTB/1000 

740  PRINT  #1,"  FC  (PSI)  -",UC 

750  PRINT  #1,  USING  "&  ##.##"; "CONCRETE  THICKNESS  AS  CALCULATED  BY 
FAA  DESIGN,  HC  (IN)  -  ",HCD 
760  PRINT  #1, "GEAR  TYPE  :",GT$ 

765  PRINT  #1, "GROSS  LOAD  OF  DESIGN  AIRCRAFT,  P  (KIPS)  -  ", P/1000 
770  FOR  JK-1  TO  5: PRINT  #1,"":NEXT  JK 
780  UC»1091.2*EXP(3.596552E-07*ECTB) 

790  PRINT  #1,"  DETERMINATION  OF  CONCRETE  AND  CTB  THICKNESS 
COMBINATIONS" 

800  PRINT  #1,"" 

805  Pl-0 

810  IF  GT$  -»SW"  OR  GT$«"sw"  THEN  Pl*84211! : J-l:GOTO  840  ELSE  GOTO 
820 
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820  IF  GT$  «"DW"  OR  GT$»"dv"  THEN  Pl*190000! : J«2 :GOTO  840  ELSE 
GOTO  830 

830  IF  GT$  »"DT"  OR  GT$-"dt"  THEN  Pl-305263 !  :  J-3  :GOTO  840  ELSE 
GOTO  660 

840  SNB=P/P1*(A(J,7)+A(J,1)*L0G(ESUB)+1/HCDA2*(A(J,3)*L+A(J,4)*K) 
+LOG(L)*(A(J,5)*HCD+A(J,2)*LOG(ESUB) )+(A(J,6) *(1/HCDA3) ) ) 

850  GOTO  950 
860  FOR  1=1  TO  3 
870  FOR  J=1  TO  7 
880  READ  A(I, J) 

890  NEXT  J 
900  NEXT  I 

910  DATA  79.87059, -20.83124 , 1043 . 40315, - 
28.99054,1.11244,33382.69811,-54.80360 
920  DATA  70.87699,-13.20461,1886.23113,-35.41080,0.91015,0.0,- 
307.15867 

930  DATA  -104.54994,41.08126,1754.38235,31.44690,-1.86926,- 
74211.45618,-601.92659 
940  RETURN 
950  GOSUB  1400 
960  GOSUB  1630 

970  COUNT  -  0:HL-0:HO0:HCTB-0 
980  HL  -  HL+10 
990  HC-HC+8 
1000  HCTB-HCTB+4 

1005  PRINT  #1,"  THICKNESS  (IN)  STRESS  (PSI) 

STRESS  ALLOWABLE" 

1010  PRINT  #1, "CONC.  CTB  PREDICTED  DESIGN  RATIO 

PASSES" 

1020  PRINT  #1,"" 

1030  ER1-4 000000 J/ECTB 
1040  ER2-ECTB/ESUB 
1050  COUNT-COUNT+1 

1060  Ll-( (4000000! *HCA3)/ (11. 73*K) ) A .  25 
1070  L2-((ECTB*HCTBA3)/(11.73*K))A.25 

1080  SCB1=B(J,13)+B(J,1)*L0G(L1)/L1+HC*L0G(L1)*(B(J,2)+B(J,5)*L2) 
+1/HCA2* (B(J, 3) *L1+B(J, 6) *ER1*L1+B(J, 9) *L2+B(J, 11) *K* (L1/L2) + 
B( J, 8) *L2*L1) 

1090  BB=1/HCA . 33*L0G(L1/HC) * (HCTB/HC) A2*ER1A-. 25 

1100  SCB2-BB* ( B ( J , 4 ) +B ( J , 7 ) *ER2+B (J, 10) *L2+B(J, 12) * (ER1/ESUB) A . 2 5 ) 
1110  SCB=P/P1* (SCB1+SCB2) 

1120  IF  COUNT  >1000  THEN  GOTO  1380 

1130  IF  (SCB-SNB) >2 . 5  THEN  1140  ELSE  1200 

1140  HCTB-HCTB+HL 

1150  IF  HCTB<4  GOTO  1310 

1160  HL-HL/2.5 

1170  IF  HL< . 05  THEN  HL=.05 

1180  IF  HCTB>22  THEN  HC»HC+.5:HCTB=4:HL»10 

1190  GOTO  1050 

1200  IF  (SNB-SCB) >2 . 5  THEN  1250  ELSE  1210 
1210  GOSUB  1520 

1220  PRINT  #1, USING  "##.#  ##.#  ####  #### 

#.###  ##.##AAAA";HC,HCTB,SCB,SNB,SR,N 
1230  IF  HC>26  GOTO  1310 
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1240  JhC-riC+ .  5 : HCTB-4 : HL-10 : GOTO  1050 

1250  HCTB-HCTB-HL 

1260  IF  HCTB<4  GOTO  1310 

1270  KL-HL/2.5 

1280  IF  HLc.05  THEN  HL-.05 

1290  IF  HCTB>22  THEN  HC-HC+.5:HCTB-4:HL=10 

1300  GOTO  1050 

1310  PRINT  #1,"" SPRINT  #1,"» 

1320  PRINT  #1, "FOR  THE  LAST  ITERATION:" 

1330  PRINT  #1, USING  "i  ##  4*;"HCTB  IS  <  4  IN.  -  USE  HCTB*  4  IN. 
AND  HO",  HC,  "IN. " 

1340  IF  C$»"S"  OR  C$-"s"  THEN  1360  ELSE  1350 
1350  PRINT  #1, CHR$ (12) 

1360  CLOSE  #1 
1370  END 

1380  PRINT  "NO  CONVERGENCE:  CHANGE  MODULUS  OF  CTB" 

1390  END 

1400  FOR  1*1  TO  3 
1410  FOR  KK-1  TO  13 
1420  READ  B(I,KK) 

1430  NEXT  KK 
1440  NEXT  I 

1450  DATA  3152 . 60047 ,0.74036, 1206.42011,40.11704 , 2 . 891281e- 
03,3.80009,0.02337 

1460  DATA  0.0, -925. 25867, 0.45975, -5. 58656, -44. 92398, -30/. 75101 
1470  DATA 

3903.81004,1.06777,1839.96268,79.37649,0.0,6.99535,0.03171,0. 

0 

1480  DATA  -1421.66496,0.73535,-10.42594,-195.49499,-388.66280 
1490  DATA  543.62830,0.0,1304.04836,66.63757,0.0,5.63003,0.03013,- 
5.30015 

1500  DATA  -820.59149,0.36554,-7.26875,-139.49406,-41.35524 
1510  RETURN 

1520  SB1- (ER1/ESUB) *.25*(C(J,5) *HC*LOG ( LI ) +C ( J , 6) *ER2+C ( J , 10 ) *ER1+ 
C(J,12)) 

1530  SB2*K/HCA2*(C(J,1) *HCTBA2+C(J, 4) *L1/L2) 

1540  SB3*L0G(L1) *HC* (C(J, 3) *L1/HCA2+C(J, 7) * (1/HCA . U33*L0G(L1/HC) * ( 
HB/HC) A2*ERlA-.25) +C(J, 11) *L0G(L1)/L1) 

1550  SB4*ER1*(C(J,2)+C(J ,8) *L2)+C(J,9) *L0G(L1)/L1 
1560  L0GSB*C(J, 13 ) +SB1+SB2+SB3+SB4 
1570  SB*P/P1* (EXP(LOGSB) ) 

1580  UO1091. 2*EXP(3 . 596552E-07*ECTB) 

1590  MR-UC/5 
1600  SR«SB/MR 

1610  N«10A ( ( . 972-SR) /8 . 250001E-02) 

1620  IF  SR<. 5  THEN  GOTO  1220  ELSE  GOTO  1220 
1630  FOR  1*1  TO  3 
1640  FOR  N-l  TO  13 
1650  READ  C(I,N) 

1660  NEXT  N 
1670  NEXT  I 

1680  DATA  -1.54962E-04, -0.09001, 0.32334, -0.02994, -0.0185, - 
6.82371E-03 


1690  DATA  -7. 93621E-03, 2. 778992E-04, 22. 59014, 0.47621, -0.15486,- 
22.51655,1.57794 

1700  DATA  -1. 14871E-04 , -0 . 08893 ,0.3387,-0.0235,-0. 01829 , -7 . 46144E- 
03 

1710  DATA  -7. 78525E-03, 5. 404066E-04, 20. 06895, 0.46020, -0.12404,- 
22.89465,1.92671 

1720  DATA  -9. 69983E-05, -0.09306, 0.36498, -0.02921, -9. 29258E-03,- 
8.35922E-03 

1730  DATA  -7. 471E-03, 7. 702747E-04, 17. 86815, 0.45719, -0.08699, - 
23.87935,1.46192 
1740  RETURN 
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CTBEVL  PROGRAM 


10  KEY  OFF: CLS 
20  LOCATE  10,5 

30  PRINT  "THIS  PROGRAM  PROVIDES  A  LIST  OF  THICKNESS  COMBINATIONS 
OF  CONCRETE" 

40  PRINT  "  SURFACE  AND  CEMENT  TREATED  BASE.  THE  SELECTED 
THICKNESS  COMBINATIONS" 

50  PRINT  "  ARE  BASED  ON  ACHIEVING  AN  EQUIVALENT  MAXIMUM  TENSILE 
STRESS  AT  THE" 

60  PRINT  "  BOTTOM  OF  THE  CONCRETE  SURFACE  LAYER.  STRESS  RATIOS 
IN  THE  BASE  ARE" 

70  PRINT  "  PROVIDED  FOR  EACH  COMBINATIONS.  A  STRESS  RATIO  OF 
0.5  OR  LESS  IS" 

80  PRINT  "  RECOMMENDED." 

90  LOCATE  24,25 

100  PRINT  "PRESS  ANY  KEY  TO  CONTINUE" 

110  A$»INKEY$ 

120  IF  A$-""  THEN  GOTO  110  ELSE  130 

130  CLS 

140  LOCATE  5 

150  PRINT  "INPUT  REQUIREMENTS  INCLUDE  SUBGRADE  PROPERTY,  THICKNESS 
OF  CONCRETE" 

160  PRINT  "SURFACE  LAYER  AS  DETERMINED  BY  FAA  DESIGN  PROCEDURE, 
LOAD  CONDITION," 

170  PRINT  "AND  THE  CEMENT  TREATED  BASE  MATERIAL  PROPERTY.  THE 
STRESSES  USED" 

180  PRINT  "TO  ACHIEVE  THE  EQUIVALENCY  BETWEEN  THE  TWO  PAVEMENT 
SYSTEM  ARE" 

190  PRINT  "COMPUTED  USING  LINEAR  REGRESSION  EQUATIONS.  THE  VALID 
RANGE  OF" 

200  PRINT  "VARIABLES  IS  LISTED  BELOW.  NO  ERROR  MESSAGES  WILL  BE 
•PRINTED  IF" 

210  PRINT  "THE  SELECTED  VARIABLE  VALUE  IS  OUTSIDE  THE  VALID  RANGE 
OF  THE  EQUATIONS." 

220  LOCATE  13 

230  PRINT  "  SUBGRAGE  RESILIENT  MODULUS  =  5  TO  30  KSI" 

240  PRINT  "  MODULUS  OF  SUBGRADE  REACTION  =  60  TO  240  PCI" 

250  PRINT  "  THICKNESS  OF  CONCRETE  LAYER  =  8  TO  26  INCHES" 

260  PRINT  "  THICKNESS  OF  CTB  =  4  TO  22  INCHES" 

270  PRINT  "  LOAD  CONDITION  “SINGLE  WHEEL  GEAR,  DUAL  GEAR," 

280  PRINT  "  AND  DUAL  TANDEM  GEAR" 

290  PRINT  "  CTB  RESILIENT  MODULUS  =  250  TO  2,500  KSI" 

300  PRINT  "  CTB  COMPRESSIVE  STRENGTH  =  1,200  TO  2,700  PS I" 

310  LOCATE  24,25 

320  PRINT  "PRESS  ANY  KEY  TO  CONTINUE" 

330  A$-INKEY$ 

340  IF  A$“""  THEN  GOTO  330  ELSE  350 
350  CLS 

360  OPTION  BASE  lsDIM  A(3,7),  B(3,13),  C(3,13) 

370  GOSUB  860 
380  LOCATE  10 


C7 


390  INPUT"DO  YOU  WANT  TO  SEND  OUTPUT  TO  SCREEN,  PRINTER,  OR  FILE 
(S,  P,OR  F) " ?C$ 

400  IF  C$«"S"  OR  C$**"sn  THEN  430 

410  IF  C$«"f"  OR  C$*»"F"  THEN  GOTO  440 

420  IF  C$«"p"  OR  C$-»P"  THEN  GOTO  460  ELSE  390 

430  OPEN  "SCRNs"  FOR  OUTPUT  AS  #l:GOTO  470 

440  INPUT "NAME  OF  YOUR  OUTPUT  FILE" ;STT$ 

450  OPEN  STT$  FOR  OUTPUT  AS  #l:GOTO  470 
460  OPEN  "LPT1"  FOR  OUTPUT  AS  #1 
470  CLS 
480  LOCATE  10 

490  INPUT"WHICH  SUBGRADE  PROPERTY  DO  YOU  WANT  TO  USE  (K  OR  E)"?A$ 
500  IF  A$«"K"  OR  A$="k"  THEN  GOTO  540 

510  IF  A$»"E"  OR  A$="e."  THEN  GOTO  520  ELSE  GOTO  490 

520  INPUT "ENTER  SUBGRADE  RESILIENT  MODULUS,  E 
(KSI) " ;ESUB:ESUB=ESUB*1000 

530  K*10A(((LOG(ESUB)/LOG(10))-1.415)/1.284) s  GOTO  560 
540  INPUT "ENTER  MODULUS  OF  SUBGRADE  REACTION,  K  (PCI)";K 
550  ESUB«10A (1.284* (LOG(K)/LOG(10) )+l. 415) 

560  INPUT"WHICH  CTB  MATERIAL  PROPERTY  DO  YOU  WANT  TO  USE  (FC  OR 
E)"?B$ 

570  IF  B$=>"FC"  OR  B$»"fc"  THEN  GOTO  610 

580  IF  B$-"E"  OR  B$»"e"  THEN  GOTO  590  ELSE  GOTO  560 

590  INPUT " ENTER  MODULUS  OF  CTB,  E  (KSI) " ;ECTB: ECTB=ECTB* 1000 

600  UC»1091.2*EXP(3.596552E-07*ECTB) :GOTO  630 

610  INPUT" ENTER  COMPRESSIVE  STRENGTH  OF  CTB  (PSI)";UC 

620  ECTB»(LOG(UC)-LOG(1091.2))/3.596552E-07 

630  INPUT"ENTER  CONCRETE  THICKNESS,  HC  (IN.),  AS  DETERMINED  FROM 
FAA  DESIGN" ?HCD 
640  ER-4 000000 1/ESUB 

650  Ir3 (4000000!  *  (HCDA3) /  (11 . 73 *K) )  A .  25 
660  INPUT "ENTER  GEAR  TYPE  (SW, DT, DW) " ;GT$ 

665  INPUT "ENTER  AIRCRAFT  GROSS  LOAD,  P  (KIPS) " ;P: P=P*1000 
670  CLS 

680  IF  C$-"p"  OR  C$="P"  OR  C$»"F"  OR  C$«"f"  THEN  GOTO  690  ELSE 
GOTO  710 
690  LOCATE  12,15 

700  PRINT  "PLEASE  WAIT  FOR  OUTPUT  FILE  OR  PRINTER" 

710  PRINT  #1, "SUBGRAGE  PROPERTIES  -  ESUB  (KSI)  =  »,ESUB/1000 

720  PRINT  #1,"  K  (PCI)  «  ",K 

730  PRINT  #1, "BASE  PROPERTIES  -  ECTB  (KSI)  =",ECTB/1000 

740  PRINT  #1,"  FC  (PSI)  »",UC 

750  PRINT  #1,  USING  "&  ##.##"; "CONCRETE  THICKNESS  AS  DETERMINED 
FROM  FAA  DESIGN,  HC  (IN)  =  ",HCD 
760  PRINT  #1, "DESIGN  GEAR  TYPE  =",GT$ 

765  PRINT  #1, "GROSS  LOAD  OF  DESIGN  AIRCRAFT,  P  (KIPS)  =  ",  P/1000 
770  FOR  JK=1  TO  5 : PRINT  #1,"":NEXT  JK 
780  UC=1091.2*EXP(3.596552E-07*ECTB) 

790  PRINT  #1,"  DETERMINATION  OF  CONCRETE  THICKNESS  FOR  A  GIVEN 
BASE  THICKNESS" 

800  PRINT  #1, ""SPRINT  #1,"" SPRINT  #1,"" 

810  IF  GT$  ="SW"  OR  GT$="sw"  THEN  Pl=84211! S J=lsGOTO  840  ELSE  GOTO 
820 
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820  IP  GT$  *»"DW"  OR  GT$»"dw"  THEN  Pl-190000! :J-2 : GOTO  840  ELSE 
GOTO  830 

830  IF  GT$  «"DT"  OR  GT$-"dt"  THEN  Pl-305263 l : J“3 sGOTO  840  ELSE 
GOTO  660 

840  SNB-P/P1* ( A ( J , 7 ) +A ( J , 1 ) *LOG ( ESUB) +1/HCDA 2  * ( A ( J , 3 ) *L+A ( J , 4 ) *K)  + 
LOG (L) *(A(J,5) *HCD+A (J,2) *LOG (ESUB) )  +  ( A ( J , 6 ) * ( 1/HCDA  3 ) ) ) 

850  GOTO  950 
860  FOR  1-1  TO  3 
870  FOR  J-l  TO  7 
880  READ  A(I,J) 

890  NEXT  J 
900  NEXT  I 

910  DATA  79 . 87059 , -20.83124 , 1043 . 40315 , - 
28.99054 , 1.11244,33382 . 69811,-54 .80360 
920  DATA  70.87699,-13.20461,1886.23113,-35.41080,0.91015,0.0,- 
307.15867 

930  DATA  -104.54994,41.08126,1754.38235,31.44690,-1.86926,- 
74211.45618,-601.92659 
940  RETURN 
950  GOSUB  1390 
960  GOSUB  1610 

970  COUNT  -  0 : HL»0 J HC«0  s  HCTB*0 
980  HL  =  HL+5 
990  HCSHC4*8 

1000  INPUT"ENTER  KNOWN  THICKNESS  OF  CTB" ; HCTB 

1001  CLS 

1005  PRINT  #1,M":PRINT  #1, ""SPRINT  #1,"" 

1009  PRINT  #1,"  THICKNESS  <IN)  STRESS  (PSI)  STRESS 

ALLOWABLE" 

1010  PRINT  #1,"  CONC.  CTB  PREDICTED  DESIGN  RATIO 

PASSES" 

1020  PRINT  #1,"" 

1030  ER1=4 000000 1/ECTB 
1040  ER2*ECTB/ESUB 
1050  COUNT— COUNT+ 1 

1060  Ll=( (4000000! *HCA3)/ (11. 73*K))A. 25 
1070  L2=( (ECTB*HCTBA3)/(11.73*K) ) A . 25 

1080  SCB1«B(J,13)+B(J,1)*LOG(L1)/L1+HC*LOG(L1)*(B(J,2)+B(J,5) *L2)+ 
1/HCA2* (B( J, 3) *L1+B(J, 6) *ER1*L1+B(J, 9) *L2+B(J, 11) *K* (Ll/L2)+B 
(J, 8) *L2*L1) 

1090  BB=1/HCA . 33*LOG (Ll/HC) * (HCTB/HC) A2*ER1A-. 25 

1100  SCB2=BB*(B(J,4)4B(J,7) *ER2+B( J, 10) *L2+B(J, 12) * (ER1/ESUB) A . 25) 
1110  SCB=P/P1* (SCB1+SCB2) 

1120  IF  COUNT  >1000  THEN  GOTO  1370 

1130  IF  (SCB-SNB) >1!  THEN  1140  ELSE  1200 

1140  HC=HC+HL 

1150  IF  HC<8  GOTO  1340 

1160  HL=HL/2 . 5 

1170  IF  HL< . 01  THEN  HL=.01 

1180  IF  HC  >26  GOTO  1340 

1190  GOTO  1050 

1200  IF  (SNB-SCB) >1 !  THEN  1240  ELSE  1210 
1210  GOSUB  1510 
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1220  PRINT  #1, USING  "  ##.#  ##.#  ####  ####  #.### 

##.##AAAA";HC,HCTB,SCB,SNB,SR,N 
1230  GOTO  1330 
1240  HC=HC-HL 
1250  IP  HC<8  GOTO  1340 
1260  HL-HL/2.5 
1270  IF  HL<.01  THEN  HL=.01 
1280  IF  HC>26  THEN  GOTO  1340 
1290  GOTO  1050 
1300  PRINT  #1, ""SPRINT  #1,"" 

1310  PRINT  #l,"FOR  THE  LAST  ITERATION:” 

1320  PRINT  #1, USING  "&  ##  &";"HCTB  IS  <  4  IN.  -  USE  HCTB=  4  IN. 

AND  HC=",  HC,"IN." 

1330  IF  C$="S"  OR  C$="s"  THEN  1350  ELSE  1340 
1340  PRINT  #1,CHR$ (12) 

1350  CLOSE  *1 

1360  END 

1370  PRINT  "NO  CONVERGENCE:  CHANGE  MODULUS  OF  CTB" 

1380  END 

1390  FOR  1=1  TO  3 
1400  FOR  KK»1  TO  13 
1410  READ  B (I , KK) 

1420  NEXT  KK 
1430  NEXT  I 

1440  DATA  3152. 60047, 0.74036, 1206. 42011, 40. 11704, 2. 891281e-03, 
3.80009,0.02337 

1450  DATA  0.0,-925.25867,0.45975,-5.58656,-44.92398,-307.75101 
1460  DATA  3903.81004,1.06777,1839.96268,79.37649,0.0,6.99535, 
0.03171,0 

1470  DATA  -1421.66496,0.73535,-10.42594,-195.49499,-388.66280 
1480  DATA  543.62830,0.0,1304.04836,66.63757,0.0,5.63003,0.03013,- 
5.30015 

1490  DATA  -820.59149,0.36554,-7.26875,-139.49406,-41.35524 
1500  RETURN 

1510  SB1= (ER1/ESUB) A . 25* (C ( J, 5) *HC*LOG (LI) +C ( J, 6) *ER2+C ( J, 10) *ER1+ 
C(J, 12) ) 

1520  SB2=K/HCA2* (C(J, 1) *HCTBA2+C(J, 4) *L1/L2) 

1530  SB3=L0G(L1) *HC* (C(J, 3 ) *L1/HCA2+C ( J, 7 ) * (1/HCA . 333*LOG (Ll/HC) * ( 
HB/HC) A2*ER1A— . 25) +C ( J, 11) *LOG (LI) /LI) 

1540  SB4=ER1* (C(J , 2) +C ( J, 8) *L2) +C ( J, 9) *LOG (L1)/L1 
1550  LOGSB=C ( J , 13 ) +SB1+SB2+SB3+SB4 
1560  SB=P/P1*EXP(L0GSB) 

1570  UC=1091 . 2*EXP (3 . 596552E-07*ECTB) 

1580  MR=UC/5 
1590  SR=SB/MR 

1595  N=10A ( ( . 97 2 -SR) /8 . 250001E-02 ) 

1600  IF  SR< . 5  THEN  GOTO  1220  ELSE  GOTO  1220 
1610  FOR  1=1  TO  3 
1620  FOR  N=1  TO  13 
1630  READ  C(I,N) 

1640  NEXT  N 
1650  NEXT  I 

1660  DATA  -1.54962E-04, -0.09001, 0.32334, -0.02994, -0.0185, - 
6 . 82371E-03 
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1670  DATA1  -7. 93621E-03, 2. 778992E-04 , 22 . 59014, 0.47921,-0. 15486, 
22. 51655,1. 57794 

1680  DATA  -1.14871E-04, -0.08893, 0.3387, -0.0235, -0.01829, 
-7.46144E-03 

1690  DATA  -7. 78525E-63, 5. 404066E-04, 20. 06895, 0.46020, -0.12404, 
22.89465,1.92671 

1700  DATA  -9. 69983E-05, -0.09306, 0.36498, -0.02921, -9- 29258E-03, 
8 . 35922E-03 

1710  DATA  -7. 471E-03, 7. 702747E-04, 17. 86815, 0.45719, -0.08699, - 
23.87935,1.46192 
1720  RETURN 
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